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Aims. Archival observations of infrared fine- structure lines of the young Galactic supernova remnant Cassiopeia A allow us to test 
existing models and determine the physical parameters of various regions of the fast-moving knots, which are metal-dominated clouds 
of material ejected by the supernova explosion. 

Methods. The fluxes of far-infrared [O i] and [O in] lines are extracted from previously unpublished archival ISO data. The archival 
Spitzer data are used to determine the fluxes of the O, Ne, Si, S, Ar, and Fe ion fine-structure lines originating in the fast-moving 
knots. The ratios of these line fluxes are used as plasma diagnostics. We also determine the infrared line flux ratios with respect to the 
optical [O m] 5007 A line in the knots with previously measured reddening. Additionally, we analyze several optical and near-infrared 
observations of the fast-moving knots to obtain clearer insight into the post-shock photoionized region structure. 
Results. We show that the infrared oxygen line flux predictions of all existing theoretical models are correct only to within a factor 
of a several. Comparison of the model predictions shows that to reproduce the observations it is essential to include the effects 
of the electron conductivity and dust. Detailed analysis of the diagnostic line flux ratios of various ions allows us to qualitatively 
confirm the general model of fast-moving knot emission and determine observationally for the first time the physical conditions in 
the photoionized region after the shock. We infer from the [O m] line flux ratios that the pre-shock cloud densities are higher than 
assumed in existing theoretical models and most probably correspond to several hundred particles per cm 3 . We also determine the 
Cas A luminosity in the infrared continuum and lines. We show that accounting for the charge exchange processes in the post-shock 
photoionized region allows us to reproduce most of the relevant spectral line ratios even in the frame of a single-temperature model 
of this region. We also estimate its plasma parameters, thickness, and carbon abundance. 

Key words, atomic processes - ISM: supernova remnants - infrared: ISM 



1. Introduction 

The young Galactic supernova remnant Cassiopeia A (Cas A) 
is one of only a few objects where the inner composition of 
the supernova progenitor may be studied directly, as the explo- 
sion ejecta have not yet mixed with the surrounding medium. 
In places where these ejecta are being crossed by the reverse 
shock wave, they are observable as compact knots emitting opti- 
cal forbidden lines of ions of oxygen, sulfur, argon, and other el- 
ements. Because of their high density, these ejecta knots are not 
strongly decelerated by the reverse shock and, as a result, have 
high proper motions that was the reaso n to name them as the fast- 
movin g kn ots (FMKs). Studies o f iPeimbert & van den Bergh 
d!97lh and IChevalier & Kirshnen 1 d 19781 1 19791) firmly estab- 
lished that the FMKs are an observational manifestation of the 
supernova ejecta. 

In FMKs, we therefore deal with a very unusual situation in 
astrophysics, where oxygen comprises 80%-90% of all atoms in 
the plasma, other components being mostly Ne, Si, S, Ar, and 
Fe. The hydrogen and helium abundances are, in contrast, un- 
derstood to be negligibly low. Optical observations have shown 
that some knot s consist even of almost pur e oxygen (so-called 
[O m] filament, IChevalier & Kirshneii d 19791) ') . 

Theoretical interpretation of the optical observations has 
proven that the plasma heated up to temperatures of several mil- 
lion degrees is rapidly cooling because of its high density and the 
anomalous chemical composition of the gas. At the same time, 
the ionization degree lags behind the rapid cooling rate. This re- 



sults in a unique situation where atoms with strongly different 
ionization degrees coexist at low temperatures from 10 4 K to 
several hundreds Kelvin. Because of this, as we now show, the 
lines of highly-charged ions, such as [Six] and [Arv], originate 
at temperatures of about 10 4 K and loweiQ. 

Traditionally, theoretical models of fast-moving knot emis- 
sion were compared only with optical line observations. Since 
there is only a very limited number of diagnostic line ratios in 
the optical spectrum, the FMK theoretical models are difficult to 
both construct and compare with optical observa tions . The m ost 
well-known models are t hose constructed bv | ltohl d 1981 alibi) . 
iBorkowski & Shulll (fl99(il), andlSutherland & Dopital dl995l) 

iDocenko & Sunvaevl d2008l) illustrated one way of constrain- 
ing the theoretical models using prospective observations of op- 
tical and near-infrared recombination lines of highly-charged 
oxygen ions. 

In this paper, we demonstrate how existing archival in- 
frared observations can discriminate between existing models 
and place strong constraints on the construction of future ones. 
We compare the model predictions of the ratios of far-infrared 
(FIR) fine-structure oxygen line fluxes to the flux in the opti- 
cal [O ra] 5007 A line with their observational values, and show 



1 The [Six] l ine has also been detected in spectra of several active 
galactic nuclei (Thompson. 119961 : iRiffel et all 12006). where it is ap- 
parently formed in the photoionized region. Even so, the models pre- 
dict line-forming region temperatures of around 10 5 K (Ferguson et al., 
Il997h . 
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that all of the predictions are incorrect by a factor of several for 
one or several flux ratios. It should be remembered, however, 
that the theoretical models predict the deviations of ionic abun- 
dances from their collisional equilibrium values by several or- 
ders of magnitude. Therefore, inconsistencies of the observed 
order suggest that some corrections should be applied to the ex- 
isting models, while the general picture described in the models 
is entirely correct. 

The infrared fine-structure line flux ratios of other element 
(in addition to oxygen) ions are currently the most reliable tool 
available for direct studies of the physical parameters in the 
FMKs. We attempted to estimate these parameters using avail- 
able observational data. As a result, parameters of the post-shock 
photoionized region were estimated observationally for the first 
time, confirming qualitatively the predictions of the theoretical 
models, but again highlighting some quantitative differences. 

The FIR lines will provide even more important probe in the 
future, as several far-infrared observatories, such as the Herschel 
Space Observatory and SOFIA, begin their operation. These 
telescopes will provide far higher quality FIR data ofhigher sen- 
sitivity and angular resolution, allowing us to constrain the FMK 
models more tightly. Results of this present study are also impor- 
tant to understand the small-scale structure of other oxygen-rich 
supernova remnants, such as Puppis A, N132D, and G292+1.8. 

The paper structure is the following. In the next section, 
we describe briefly the existing theoretical models of the FMKs 
and estimate oxygen infrared line fluxes in one of these models. 
Section [3] describes the archival infrared observations of Cas A 
and their analysis. In Sect. [4] we estimate the physical conditions 
and abundances in different regions of FMKs from the observed 
line flux ratios. We also discuss these results and compare them 
with values predicted by the models. In Sect. [5] we analyze spec- 
tral line ratios to obtain a consistent model of the post-shock 
photoionized region. Section [6] is devoted to the predictions of 
the mid- and far-infrared recombination line intensities. Finally, 
in Sect. [7] we present our conclusions. 



2. Theoretical models of fast-moving knots 

The first detailed model of fast-moving knot emission, describ- 
ing it as originating from the passage of a shock th rough the pure 
oxygen medium, was constructed by lltohl (11981 alibi) . The predic- 
tions of his model H were compared with measurements of the 
[O m] filament optical spectrum in the northern part of the Cas A 
supernova remnant (SNR) shell, and were found to reproduce all 
four measured oxygen [Oi], [On], and [Oin] optical line flux 
ratios to within a factor of two. The infrared line ratios to the 
[O m] 5007 A lin^| for this theoretical model (denoted I-H) and 
other models described below are given in Table Q] 

In this and other theoretical models, the rapidly cooling re- 
gion just after the shock front produces high ionizing radiation 
flux that in its turn produces two photoionized regions (PIRs), 
bef ore a nd after the shock front (see Fig.[TJ. 

lltohl (Qj>86) found that the model neutral oxygen optical line 
intensities originating in the PIR after the shock front are much 
too high compared to the observed spectra of the oxygen-rich 
supernova remnant Cas A and Puppis A ejecta. He suggested 

2 The 5007 A line is three times brighter than the other [O in] doublet 
component at 4959 A (e.g jQsterbrock & Feriancj i2006l) 1. Note that in 
some papers the line intensities are compared to the sum of the doublet 
components. 
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Fig. 1. Schematic representation of the FMK temperature struc- 
ture, induced by its interaction with the reverse shock. Negative 
column densities corresponding to the pre-shock regions are in 
linear scale; positive (post-shock) column densities are on a log- 
arithmic scale. The post-shock photoionized region (PIR) pa- 
ramet ers are given as derived in S ect.[5](solid line) and according 
to the lBorkowski & Shullld 19901) models (dashed line). Shock is 
moving to the left. Because of the FMK high density, the reverse 
shock in it has slowed down to about 200 km/s. 



that emission from this region is damped because the region it- 
self is truncated due to some hydrodynamical phenomena. This 
truncation affects essentially only neutral oxygen lines, as all 
the other ions have already recombined in the dense post-shock 
plasma before reaching the photoionized region. This possibly 
ensures that the [0 1] lines are the least reliable for direct compar- 
ison with model predictions, since all existing models are one- 
dimensional. 

Another group of theo retical FMK models was proposed by 
iBorko wski & Shul (I1990I) . In contrast to the Itoh models, some 
of these accounted for electron conductivity, which changed 
the line fluxes consider ably. One more difference is that the 
IBorkowski & Shulll (1 19901) models do not account for the emis- 
sion from the PIR before the shock. The models that most accu- 
rately describe the Cas A FMK optical spectra are F and DC (re- 
spectively, denoted as BS-F and BS-DC in Table Q] and below). 
We note that only the BS-DC model takes the electron conduc- 
tivity into account. 

Both ||toh| (Il981allbl) and IBorkowski & Shulll (1 19901) assume 
that the emission originates in t he pur e oxygen plasma. However, 
it was shown by Dopit a~et~aTl ( 1 19841) that inclusion of other el- 
ements in the model significantly changes the plasma thermal 
structure and emission, especially in the cold PIRs, where a 
number of fine-structure transitions in ions of other elements 
can effectively cool plasma to temperatures of about one hun- 
dred Kelvin. Following this line o f reasoning, another group o f 
the FMK models was published bv lSutherland & Dopital (fl995). 
These models assume that the emission is produced by the in- 
teraction of the dense cloud with the external shock, entering 
the cloud and propagating through it. The expected far-infrared 
line intensities are not published and we estimate them in the 
Sects. 12. 1 1 and [2721 for the model with a 200 km/s shock speed 
and pre-shock ion number density of 100 ctrT 3 (denoted as SD- 
200 in Tableland below), which most accurately reproduces the 
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Table 1. Far-infrared line flux ratios to the dereddened 5007 A line as derived from different theoretical models and observations. 
See text for details. 





A, yum 






///(5007) 




/pir 






I-H 


BS-F 


BS-DC 


SD-200* 


Observed** 


SD-200* 


[Oiv] 


25.91 


0.0051 


5.0 


0.53 


8.1 


0.20 


0.995 


[Out] 


51.81 


0.031 


1.34 


0.72 


0.67 


0.25 


0.999 


[Oi] 


63.19 


1.30 


170 


12 




0.07 




[Out] 


88.36 


0.028 


0.22 


0.11 


0.55 


0.10 


0.999 


[Oi] 


145.5 


0.038 


7.8 


0.75 




<0.0024 





* Derived in Sects. l2.1l and l2.2l 
** Obtained in Sects. j37JJ and E2J] 



Notes. The last column lists contribution to the line flux in the SD-200 model arising from the photoionized region (PIR) before the shock. 
Derivation of the observed line flux ratio values is described in Sect. [3] The model-predicted [Oi] lines are sensitive to the pre-shock ion number 
density in the model, see Sect. 14.1.21 



optical spectra of Cas A fas t-moving knots, as demonstrated by 
ISutherland & Dopital dl995l) . 

Compared to the models of lltohl d!981albl) and 
iBorkowski & Shulll (Il990l) . the SD-200 model does not in- 
clude the PIR after the shock front, because it is constructed 
to describe the FMK optical spectrum and the temperature of 
this post-shock PIR was estimated to be too low to contribute 
significantly to the optical line emission. 

All described models nevertheless share many similar fea- 
tures (see Fig. [TJ. In all of them, the plasma after the shock 
front passage rapidly cools and, at temperatures of (5 - 50) x 
10 3 K, emits the high-ioniz ation lines observable in the visi- 
ble and near-infrared spectra dChevalier & Kirshnerl 119781 119791; 
iHurford & FesenLll996tlGerardv & FesenLl200ll) . The thickness 
of this emitting layer is extremely small (about 10 10 cm for the 
pre-shock ion density of 100 cm -3 ), but because of the high elec- 
tron density in the cooling material at T < 10 5 K (more than 
10 5 cm" 3 ), the emission measure is enough to produce or con- 
tribute to bright emission lines of highly-charged ionic species, 
such as [O m] lines near 5000 A. 

In all of these models, the visible line emission of weakly- 
ionized ionic species (e.g., [Oi] and [On]) originates in pho- 
toionized regions before and after the shock front. 

Although all of the theoretical models describe optical spec- 
tra relatively well, with discrepancies in relative intensities that 
do not exceed factors of about two, predictions for the infrared 
lines differ by several orders of magnitude (see Table Q}. It is 
therefore clear that the models have a lack of constraining power 
and more diagnostic information in the form of various line flux 
ratios is needed to determine the true structure of the fast- moving 
knots. Part of this information may be obtained from the far- 
infrared line archival observations, and we present this analysis 
in Sect. [3]below. 

2.1. Infrared lines from the SD-200 model 

Because of a lack of published values, we had to make our own 
estimates of the far-infrared line fluxes and their ratios with re- 
spect to the [O m] 5007 A line in the SD-200 model. Fortunately, 
all the data necessary for this calculation - the ionic abundances 
of all oxygen ions, the temperature and density structure of 
the post-shock region, and the temperature structure of the pre- 
shock photoionization front - are knowifl 



3 We performed our own computations of the post-shock plasma re- 
combination and discovered that our oxygen ion distribution over ion- 
ization stages is rather similar to that presented in the lower panel of 



We note that the FMK structure calculations in the SD-200 
model are based on the plasma composition that differs slightly 
from the composition inferred from the X-ray and optical obser- 
vations. Specifically, it is dominated by O (63% by number) and 
Ne (29%), and contains minor amounts of C, Mg and Si. In con- 
trast, the observations imply that the majority of the currently 
bright X-ray and optical plasma is oxygen (about 90%), other 
abundant elements being not only Ne, but also Si and S. 

In Fig. [2] we demonstrate how the ionization state distribu- 
tion in the SD-200 model st rongly differs from the collisional 
ionization equilibrium (CIE, Maz zotta et al.l dl998l) ) on exam- 
ples of the 2+ and 3+ ions producing bright FIR lines. It is 
seen that in the rapidly cooling plasma these ions are abundant 
at significantly lower temperatures (down to 300 K and below) 
than in the CIE because of very rapid plasma cooling by means 
of line emission. 

In the same figure, we present the oxygen ion emission mea- 
sure distribution over temperature dEo/d(lnT e ) in the SD-200 
model, defined by Eq. Q. This distribution together with the 
line emissivity dependence on temperature allows one to calcu- 
late the relative contributions to the total line flux from different 
temperature intervals (see Sect. l2.2l >. 

2.2. Line flux computation 

The emissivity, cm 3 /s, of a spectral line corresponding to transi- 
tion from level u to level I is defined as 

s(u,l) = n u A u i/n e n i7 (1) 

where « e , «i, and n u are the number densities of electrons, the 
corresponding ionic species, and ions in the state u, respectively, 
and A u i is the spontaneous u — > I transition rate. 

For oxygen ions at temperatures above 1000 K, the 
line emis sivities g were obtained from the Chianti atomic 
database dDere et all Il997t lLandi et all |2006). For lower tem- 
peratures and other ions, we computed emissivities using a stan- 
dard approach, i.e., calculating the level populations by solv- 
ing a system of linear equations describing transitions between 
the lowest ionic levels using electronic excitation effective col- 
lision strength data from the literature and spontaneous transi- 

Fig. 3 of the Sutherland & DoDi3 i 19951) . but only if the ion spectro- 
scopic symbols in th at figure are increased by u nity. Therefore, we 
assume that Fig. 3 of Sutherland & DoDita l l 19951) has a misprint and 
the ion spectroscopic symbols should be read, e.g., "O VII" instead of 
"O VI", "O VI" instead of "O V", etc. From this, it follows that the 
SD-200 model values for the neutral oxygen abundances are not known 
and we do not provide estimates of the FIR [O i] line intensities. 
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Fig. 2. Comparison of ionic abundances of oxygen ions 2+ and 
3+ in the collisional ionization equilibrium (thin curves) and 
the cooling post-shock region of the SD-200 model (thick curves 
of the same color). The differential oxygen emission measure per 
logarithmic temperature interval dEo/dlnT e is indicated by the 
red solid line. 



tion rates from both the Ralchen ko et al.l (|2007) database^ and 
MCHF/MCDHF collection^, and then applying Eq. © to derive 
the line emissivities. 

The electronic collisional e xcitation effective c ollision 
strength values were t a ken from Dufton & Kingston! (I1991I), 
Blum & P radhanl (1 1992b. IZhang et all (1 19941). iButler & Zeippenl 
(119941). iLennon & Burkd d 19941). iGalavis et all (Il995h . 
Pelan & Berringtonl (1 19951). iBerrington et al l (11998b, 
Taval & Gupta! (11999). Griffin & Badnelf (2000), GriffinetdJ 
(20011). iTaval (120061) . and iRamsbottom et al.l (120071) and 
extrapolated to lower temperatures by a constant, if needed. 

Line fluxes I(u, /), erg/cm 2 /s, were computed by integrating 
along the line of sight 



I(u, I) = hv 



AnR 2 



J" e(u, I; 



T e (r)) n e (r)rii(r)dr, 



(2) 



where hv is the photon energy, R is t he distance f rom th e ob- 
server to the emitting region (3.4 kpc. iReed et al.1 ([1995)), and 
S is the emitting region area. From this expression, it is seen 
that the line emission mostly originates in regions of the highest 
emission measure n e riodr, ionic abundance n jno, and emissiv- 
ity e (here no is the total number density of all oxygen ions). 

The integral over distance in the post-shock cooling region 
can be transformed into an integral over temperature by com- 
pleting the substitution 



dr dt , 

dr = dr e 

dt dT P 



dr 



«o,t 



. — ^shock 

dt n t 



where f s hock is the shock front speed, and n t and «o,t are the to- 
tal number densities of all ions in the plasma at a given point 
and before the shock, respectively (i.e., in the SD-200 model 
«o, t = 100 cm -3 ). It is assumed here that the photon heating is 
insignificant compared to the energy losses. 



Available online at http : //physics .nist . gov/asd3 
http : //atoms . vuse . vanderbilt . edu/ 



The parameters of these equations, i.e., the cooling function, 
and electron and ion densities as functions of temperature - were 
taken from the SD-200 model. 

The derivative df/dr e in the post-shock cooling region may 
be obtained from the energy conservation law written for one 
particle as 



A' Li ' o 

\2 nt 



+ n, 

-k K T e 



-n e A^dt - p d I — 
\n< 



where is the Boltzmann c onstant, An is the cooling function 
([Suthe rland & Dopital [T993I) . and p = (n e + n^k^T^ is the gas 
pressure. Dividing both sides by dr e , we can now obtain the 
required derivative from the expression 



3, d(« e /»t) , 3 (n e 
o^BTe— — — +- — + 1 \k B 
2 dT e 2 \ n t 



-n e A N — + — 

dr e \« t 



l 



kT e dn t 
n t dT e 



We note that in the SD-200 model, the plasma pressure is not 
constant in the post-shock cooling region (their Fig. 3) at low 
temperatures because of the contribution of the magnetic field 
pressure. 

For the purpose of qualitative analysis, we introduce the oxy- 
gen differential emission measure per logarithmic temperature 
interval 



dE 
d(lnr e ) 



«o» e dr 

dr e 



(3) 



It indicates contribution of a given logarithmic temperature in- 
terval to the total emission measure, showing where most of the 
line emission originates. Using this notation, we can express the 
line flux given by Eq. (T3J as 



/(«, = hv 



4nR 2 



J s(u, I; 



dE 
d lnT e 



n\ dr e 
n T e 



Contributions from different temperature intervals to the emis- 
sion of the post-shock cooling region in the oxygen lines are 
illustrated in Fig. [3] where we plot the normalized line fluxes per 
logarithmic temperature interval. It is seen that the post-shock 
contribution to the visible [O m] line emission originates at tem- 
peratures of (5 - 20) x 10 3 K, whereas the infrared lines are emit- 
ted mostly at lower temperatures. A far broader range of tem- 
peratures also contributes to the line emission because of lower 
excitation energies. 

In our analysis, we artificially stop integrating Eq. (f2]i in 
the post-shock region when the plasma temperature drops be- 
low 300 K, since the SD-200 model does not contain data for 
lower temperatures. This should not influence the total FIR line 
intensities much, since the model predicts that they are mainly 
emitted in the pre-shock region. 

In the pre-shock PIR, the ion density is constant and the tem- 
perature dependence on time for the SD-200 model is known ex- 
plicitly. Therefore, Eq. (0 was applied directly to compute the 
line fluxes from this region. The resulting line flux ratios are 
given in Table [TJ 

3. Archival observations of the FIR lines 

Th e first search for the fO il 63.19 pm line was performed 
by iDinerstein et al.l d 1987b . who observed the northern part of 
the Cas A using the NASA Kuiper Airborne Observatory. This 
search resulted in a 3-<x upper limit to the reddening-corrected 
flux ratio 7(63. 19 pm)/ 1(5007 A) of about 0.3, which tightly con- 
strained some theoretical models. Later, the far-infrared lines of 
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Fig. 3. Differential flux contribution per logarithmic temperature 
interval for the optical [Om] 5007 A line and infrared oxygen 
lines in the cooling post-shock region according to the SD-200 
model. Curves are normalized so that the area enclosed under 
each of them equals one. Low-temperature (below 4000 K) con- 
tribution to the optical [O m] line is due to recombination. 



O, Ne, Si, S, and Ar ions were detec ted by the Infrared Space 
Observatory (ISO. lArendt et al.l[l99 9). We describe our analysis 
of some of these archival observations in Sect. I3.ll The Spitzer 
Space Telesco pe recently spectrally mapped the Cas A super- 
nova remnant dEnnis etall l2006t iRho et al.L I2008I) . In Sect. [3721 
we determine the line ratioJJofO, Ne, Si, S, Ar, and Fe lines to 
the [Ora] 5007 Aline. 

To determine these ratios, we compare the infrared observa- 
tions of ISO and Spitzer with the optical images of the Hubble 
Space Telescope. The most accurate way of comparing these ob- 
servations would be an analysis of the properties of individual 
knots but, unfortunately, the ISO data are of insufficient angular 
resolution (about 40")- The angular resolution of the Spitzer data 
is superior (2"-8"), but also too coarse to isolate the contribu- 
tions of individual knots distinguishable within the optical maps 
(072-1", corresponding to (3 - 15) x 10 16 cm at 3.4 kpc). 

As the optical obser vations show, the bright FMK lif etime is 
of the order of 30 years ( Kam per & van den BerehLI 19761) . which 
is determined by the time needed by the shock to cross the cloud. 
After that, the knots fade, but others gradually appear in the maps 
of the Cas A. Therefore, one may directly compare either obser- 
vations performed with a short time difference between them (as 
in the case of Hubble Space Telescope observations compared 
with Spitzer data), or averaged over a large area (as in the case 
of the ISO data). 

3.1. ISO observations 

The Infrared Space ObservatorjQ Long Wavelength 
Spectrometer (ISO LWS) observations of the Cassiopeia A 
supernova remnant were performed in 1996 and 1997. Here we 



6 Here and everywhere below the "line ratio" denotes the ratio of the 
line fluxes in energy units. 

7 Results in this section are based on observations with ISO, an ESA 
project with instruments funded by ESA Member States (especially 
the PI countries: France, Germany, the Netherlands, and the United 
Kingdom) and with the participation of ISAS and NASA. 




Fig. 4. ISO LWS apertures (green circles) overlaid on the Hubble 
Space Telescope 2 004 ACS image of Cas A in F475W filter 
dFesen et al.L 12006). North is upwards and East is to the left. 



Table 2. Summary of the ISO LWS archival observations of the 
Cassiopeia A. 



ID 


Date 


Exposure 


R.A. 


Dec. 


Cas A #1 


1996 Jun 24 


1054 s 


350.8656 


+58.8130 


Cas A #2 


1996 Jun24 


1052 s 


350.8614 


+58.8361 


Cas A #3 


1996 Jun 24 


1054 s 


350.8279 


+58.7919 


Cas A #4 


1996 Jun 24 


1054 s 


350.9748 


+58.8400 


Cas A #5 


1997 Jun 09 


1612 s 


350.9097 


+58.8228 


Cas A #6 


1997 Jun 09 


1612 s 


350.8963 


+58.7939 


Cas A #7 


1997 Jun 02 


1614 s 


350.8186 


+58.8188 



do not discuss observations at shorter wavelengths (A < 40 yum), 
since more recent data of significantly higher angular resolution 
is available from the Spitzer observatory. 

To determine the line fluxes, we used the ISO LWS cal- 
ibrated data from a highly-processed data product "uniformly 
processed LWS L01 spectra" available from the ISO data 
archival 

There are seven LWS observations in L01 mode each cov- 
ering the full instrument spectral range between 43 and 170 pm 
with medium spectral resolution (A/AA between 150 and 200). 
The regions observed by ISO are shown in Fig.|4]overlaid on the 
Hubble Space Telescope Advanced Camera for Surveys (HST 
ACS) Cas A image in filter F475W, containing [Ora] 4959 
and 5007 A lin es (Obs. ID 10286 , observations completed in 
December 2004. lFesen et all (120061) ). 

The ISO LWS observational data are summarized in Table [2] 
and the spectral cuts containing the discussed lines are presented 
in Fig. |5] 

To compute the far-infrared spectral line fluxes, we sub- 
tracted the background continuum flux and averaged over three 
scans of the same detector, weighting the data according to its 
uncertainties. In the case of the [Oiii] 88.36 pm line, it was pos- 
sible to estimate the systematic uncertainties by comparing the 
flux values obtained from two neighboring detectors. The flux 
differences are expected to have two causes: slightly different 
apertures of the detectors and systematic errors due to imperfect 

8 http : //www . iso . vilspa . esa . es/ida/ 
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1 45.5 |xm 



#1 #2 #3 #4 #5 #6 #7 

Fig. 5. ISO LWS spectra of Cas A at the observed locations (as shown in Fig. [4] and Table [2} are shown in the regions of oxygen 
far-infrared lines discussed in the text. Continuum emission is not subtracted. Plotted are the averages over three scans weighted 
according to their uncertainties. Two curves for the 88.36 fj.m line correspond to measurements by two different ISO detectors. The 
short-dashed line denotes the zero level; maximum intensity corresponds to 1.7 x 10~ 10 , 0.7 x 10~ 10 , 0.7 x 10~ 10 , and 0.10 x 10~ 10 
erg/cm 2 /ji/m for [O m] 51.81 /urn, [O i] 63. 19 fim, [O in] 88.36 /urn, and [O i] 145.5 fim lines, respectively. Spectral ranges are centered 
on the rest-frame wavelengths of the respective spectral lines and include a velocity range of +10 4 km/s. 



Table 3. Far-infrared line fluxes from the ISO LWS observations 
of Cas A. Fluxes from the matched regions in the HST ACS 
F475W filter are also given. The upper limits are at the 3cr level. 



Table 4. Far-infrared line flux ratios to the optical 5007 A line 
///(5007 A) corrected for the optical line reddening. 











Line 


#2 


#5 










[Ora] 51.81 yum 


0.25 


0.36 










[Oi] 63.19 yum 
[Om] 88.36,™ 


0.07 


0.10 


Line 


| #1 #2 #3 #4 #5 


#6 


#7 


0.10 


0.14 




FIR line fluxes F, 1CT 11 erg/cm 2 /s 






[Oi] 145.5 yum 


<0.0024 


<0.01 



[O m] 52 fim 
[O i] 63 yum 
[Om] 88pm 
[Oi] 145 yum 



[Om 



5.4 
<0.6 
7.5 
<0.1 



9.3 
2.6 
4.2 
<0.09 



5.4 
1.2 
2.7 
0.3 



<1.2 
<1.2 
0.5 
0.1 



3.5 
1.0 
1.9 
<0.1 



5.0 
1.0 

3.2 
<0.1 



5.2 
0.9 
3.4 
<0.1 



4959 + 5007 A line flux F, 10" 13 erg/cm 2 /s 
<0.1 25 2.0 <0.1 6.5 0.5 



2.9 



Note. ISO region positions are shown in Fig. [4] Regions #3, #6 and # 7 
have higher interstellar absorption column (Bie ging & Crutchen. 11986). 



instrumental calibration dGrv et al.Ll2003l) . Our analysis showed 
that differences between the two detector measurements do not 
exceed 10-15%, when the line flux determination accuracy is 
sufficiently high (see Fig.[5]for a visual comparison). 

The background-subtracted FIR line fluxes from the ob- 
served regions are presented in Table [3] It is seen that the flux 
at the level of up to 10% in the brightest fields is observed in 
the 88.36 yum in the region #4 situated outside the main super- 
nova remnant shell with no detected optical [Ora] emission and 
very little diffuse emission in other optical bands. We consider 
the flux from this region as originating in the foreground and/or 
background of Cassiopeia A and in our analysis subtract it from 
the fluxes of the other regions. 



Note. The 5007 A line contribution of 75% to the total FMK emission in 
the F475W filter is assumed. Optical line attenuation due to reddening 
of 200 times is assumed in both the ISO LWS regions, see text. FIR line 
attenuation is neglected. 



In the same table, we indicate the total fast-moving knot op- 
tical flux for the matched regions in the F475W filter of the dis- 
cussed HST ACS data. The diffuse emission in this filter is as- 
sumed to come only from the [O m] 4959 and 5007 A lines. 

From these data, we calculated the line ratios to the [O in] 
5007 A line, given in Table |4] Assuming that the FIR line emis- 
sion observed by ISO originates in the same spatial regions as the 
optical line emission, these flux ratios are characteristic of the 
knot emission and may be used for tests of the theoretical mod- 
els. If significant large-scale diffuse emission in the FIR lines is 
present, as observed in the [S m] and [O iv] lines by Spitzer ob- 
servations (see Sect. 13.21 1. then the determined line ratios should 
be interpreted as upper limits to the values characteristic of the 
knots. 

In Table |U we provide the values corrected for the interstel- 
lar absorption of the op tical line using an average value from 
iHurford & Fesenl ([T996). This level of reddening diminishes the 
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observed visible [Om] doublet intensities by a factor of about 
200 (Ay ~ 5.0) and is characteristic o f region #2 and, possibly, 
region #5 (Bieg ing & Cratcherul986l) . 

For southern and western regions (#3, #6 and #7) the 
reddening is known to be much higher because of absorp- 
tion in the molecular clou ds with peak Ay's reaching 10-15 
(Bieging & Crutcherlll986l) corresponding to the optical line at- 
tenuation of up to factors of 10 4 and more. Therefore, we do not 
provide dereddened line ratios in these regions. In our further 
analysis in Sects. I4land l3. 31 we use only the data from the region 
#2. 

There is an eight-year time difference between ISO and HST 
observations. D uring this time, which is about 3 0% of a bright 
FMK lifetime (iKamper & van den Ber gh. 1976), some optical 
knots may have disappeared, and new ones may have appeared. 
However, averaging over a large number of optical knots con- 
tained in any ISO LWS field of view at least partially can- 
cels changes induced by the brightening or fading of individual 
knots. 

3.2. Spitzer Space Telescope observations 

The Spitzer Space Telescope spectrally mapped the 
Cassiopeia A in the infrared range using its Infrared 
Spectrograp h (IRS) producing low-reso l ution spectra from 
5 to 38 /urn (lEnnis et all 120061 : iRho etail 12008b and detected 
in this spectral band the fine-structure lines of [Oiv], [Nen], 
[Nem], [Nev], [Sin], [Sin], [Siv], [Am], [Arm], [Arv], [Fen], 
etc. 

We used basic calibrated data of these publicly available 
Spitzer observations (Program 3310) to construct t he data cube 
over t he entire remnant with the CUBISM software dSmith et all 
120071) . The background regions were de termined from the 8 yum 
MIPS images of the supernova remnant dHines et alll2004l) . 

The line flux maps were then produced using the parabolic 
approximation for the background continuum emission. This al- 
lowed the maps to be constructed reliably even for the weak lines 
in the regions of strong continuum dust emission. 

Since the Spitzer spectral maps and the HST images were ob- 
tained almost simultaneously with a time difference of only one 
month, we were able to make a direct pixel-by-pixel comparison 
of the infrared line map with the optical image. For this purpose, 
we processed the HST images by removing the stars, convolving 
with the point spread function (PSF) of a given Spitzer spectral 
map computed as described in Appendix [A] and regridding the 
HST image to match pixels to the IRS module containing each 
infrared line. These matched HST images were produced and 
compared to each of the infrared lines. 

3.2.1 . The oxygen [O iv] 25.91 pm line 

The [Oiv] 25.91 pm spectral line map of the northern Cas A 
shell region obtained from Spitzer observations is presented in 
Fig. [6] Alongside, we show the HST image in F475W filter pro- 
cessed as described above. The northern region is of primary 
importance as it has a well-measured and almost constant inter- 
stellar absorption level. 

It can be readily seen from the brightest features in Fig.[6]that 
there is a difference between the Cas A morphologies for the op- 
tical and infrared lines. The difference may have several reasons : 
variations in the interstellar reddening (Hurford & Fesen, 1996), 
emitting plasma chemical composition dChevalier & Kirshnerl 




Fig. 6. The Spitzer Space Telescope map of the northern part 
of the Cas A optical shell rich in FMKs in the [Oiv] line at 
25.91 pm. (top) compared with the HST image in F475W fil- 
ter matched in resolution and in pixels (bottom). The maps are 
centered on R.A. of 350.8614 deg and declination of 58.8361 
deg and show the region of 2.5 x 3.5 arcminut es. The red circles 
denote regions with reddening measured by iHurford & Fesenl 
dl996h . Square root scale in intensity is used to enhance weak 
features. Diffuse infrared [O iv] emission from optically-dark re- 
gions is clearly visible. 



Il979h . and FMK pre-shock density distribution over the rem- 
nant. 

In Fig.|7J we present the resulting pixel-by-pixel scatter plot 
showing interrelation between the [O iv] infrared and [O ra] opti- 
cal line fluxes. From this comparison of the maps, we derive the 
average line ratio of the optical and infrared lines, which cor- 
responds to an average reddening, FMK density, and chemical 
composition. 

In this analysis, "average" reddening was assumed to dimin- 
ish the observed optical [O m] line emission by a fac tor of 200, 
which is an average of the Hurford & Fese r3 (fl996l) measure- 
ments. As may be inferred from their paper, small-scale varia- 
tions in the interstellar reddening in this region change the opti- 
cal line attenuation within a factor of two. The observed scatter 
of the line ratio seen in Fig.|7]is of similar magnitude. 
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[O III] 5007 A dereddened flux, erg/cnr/s 



Fig. 7. Pixel-by-pixel comparison of the dereddened 5007 A line 
flux with the infrared [O iv] line map. Fluxes are given per detec- 
tor pixel, equal to (5'.'08) 2 for the IRS LL module containing the 
infrared oxygen line. Two linear fits, both with the same slope 
of 0.20, are shown with lines. Filled squares denote FMKs from 
lHurford&Fesenl(ll996l) with measured reddening. 



In the same figure, we also present data for five FMKs with 
reddening measured bv lHurford & Fesenl d 19961) . Their positions 
are shown in Fig. Q In these cases, the reddening correction 
was made according to their measurement results. Four out of 
five points lie very close to the best-fit lines, showing that a 
considerable part of the FMKs have similar excitation condi- 
tions. A remaining point (FMK 1) has relatively weaker [Oiv] 
emission, showing that the ph ysical conditions there are d iffer- 
ent. This region (Filament 1 of Baade & Minkowski d 19541) ) in- 
deed has a composition that d iffers from that of the other regions 
(IChevalier & Kir shnerl [1971 . 

In Fig. [6] infrared emission in the [Oiv] line that originates 
in the cen t ral Ca s A region with no optical nebulosity (see also 
iRho et all (120081) ) is also noticeable. In Fig. [7] this infrared- 
only emission shows up on the left as horizontal branches, cor- 
responding to an approximately constant background surface 
brightness. Two of these branches are apparent, probably corre- 
sponding to two different column densities of the 3+ ion. Pixels 
forming the lower branch are situated just outside the optical 
shell, where the optical surface brightness is already fading, but 
the infrared line emission remains present. 

The upper branch corresponds to the diffuse emission in the 
central part of the supernova remnant. Corresponding points in 
the HST image are partially outside the Fig.|7]plot range, but the 
infrared line surface brightness remains on the level of (1 - 3) x 
10~ 13 erg/cm 2 /s/pix. 

This infrared-only line emission is consistent with the exis- 
tence of extended cold photoionized regions invisible in optical 
and X-ray emission close to the center of the supernova remnant, 
possibly connected with the thermal radio absorption o bserved 
in the same regions (Kass irnet al.Lll993llDelanevll2004l) . 

We note that the [Oiv] infrared line at 25.91 yum overlaps 
with the [Fe n] line at 25.99 /im, but the contribution of the latter 
can be estimated by comparison with the [Fen] 17.94 jum line 
map. We assume that the [Fe n] 25.99/17.94 //m line ratio is con- 
stant over the remnant near the reverse shock. Then, we select a 



/([Oiv] 25.91 fim + [Fen] 25.99 [tax) 
/([Fen] 17.94 yum) 



~ 5, 



and assume it to be the upper bound to the [Fen] 25.99/17.94/im 
line ratio later used in Sect. [5] 

This bound is used to infer that the upper limit to the 
25.99 fim [Fe II] line contribution to the total flux of the observed 
line near 26 //m is about 20 — 40% in all regions of the Cas A 
used in our analysis (i.e., the Hurf ord & Fesenl d 1996b FMKs). 

3.2.2. Other infrared lines in the Spitzer spectral maps 

Similar pixel-by-pixel comparisons with the optical image were 
also performed for other bright spectral line flux maps. In 
Table |5j we provide the derived flux ratios to the optical [O in] 
5007 A line in five iHurford & Fesenl (1 19961) FMKs and the av- 
erage values of these ratios in the northern part of the Cas A 
shell. The last digit in brackets in the "average" column denote 
the observed 2-cr scatter around the average line flux ratio in 
units of the last digit. This scatter represents the true variations 
in the line flux ratios due to both chemical composition vari- 
ations (IChevalier & Kirshnerl 1 19791) and variations in physical 
excitation conditions (e.g., different pre-shock density in differ- 
ent regions of Cas A). 

From the table, it can be seen that different line flux ratios 
have different scatter around their average values. A small scatter 
in e.g., the Ar, Si, and S line ratios to the 5007 A line implies 
that the Cas A morphology in these lines is similar to the optical 
one. A larger scatter, present in e.g., Ne and Fe line ratios to 
the 5007 A line corresponds to differ ent morphologies of the 
infrared spectral maps, as also noted bv lRho et alj d2008l) . 

The obtained line ratios may also be used as tests of the theo- 
retical models and are utilized in the next section to estimate the 
average abundances and physical conditions in the bright Cas A 
knots. 

From Table |5j it is seen that the full energy flux in the in- 
frared lines detected by Spitzer exceeds one of the brightest op- 
tical line - [O m] 5007 A - and is of the same order as the total 
power in all optical lines summed up. 

3.3. Comparison of the model predictions with observations 

As the models provide predictions only for the oxygen line rela- 
tive intensities, we cannot use observations of lines of other ions 
to directly compare with the models. 

Comparing observed far-infrared line ratios with theoreti- 
cal model results (see Table [TJ, we see that the lltohl dl981allbh 
models underestimate the [Om] and [Oiv] line fluxes by fac- 
tors of 3-40, as a result of having a far lower average ion- 
ization degree after the shock passage th an in other models. 
In contrast. thelBorkowski & ShulU d 19901) model BS-F and the 
ISutherland & Dopital d 19951) model SD-200 that do not account 
for electron conductivity, both overpredict the emission in the 
[O rv] line by factors of 20 - 40. 

The remaining model (BS-DC of iBorkowski & Shulll 
dl990 )) predicts the infrared line fluxes to within a factor of sev- 
eral, except for the neutral oxygen lines, but this problem is at 
least partially resolved by increasing the pre-shock density (see 
Sect. 14.1.2b . This strongly suggests that taking into account the 
electron conductivity is essential for a model when attempting to 
reproduce the observations. 
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Table 5. Infrared line fluxes I and their ratios 7/7(5007 A) to the dereddened optical [Om] 5007 A line from the Spitzer data. The 
column "average" contains the ratio values averaged over the northern part of Cas A and their observed scatter, and other columns 
contain data for the Hur ford & Fesenl (1996) FMKs. 



Line 




Line flux /, 


10 13 erg/cm 2 /s/pixel 






Flux ratio ///(5007 A) 








FMK 1 


FMK 2 


FMK 3 


FMK 4 


FMK 5 


Average 


FMK 1 


FMK 2 


FMK 3 


FMK 4 


FMK 5 


[Am] 6.985 yimi 


3.8 


14 


17 


5.6 


11 


0.5(1) 


0.41 


0.75 


1.15 


0.46 


0.79 


[Arm] 8.991 yum 


0.46 


2.0 


1.7 


0.80 


1.5 


0.08(3) 


0.05 


0.13 


0.12 


0.07 


0.12 


[Srv] 10.51 yum 


0.24 


1.6 


1.4 


0.67 


0.91 


0.07(2) 


0.025 


0.11 


0.10 


0.06 


0.07 


[Nen] 12.81 yum 


0.31 


0.92 


0.89 


1.0 


0.30 


0.10(8) 


0.04 


0.07 


0.07 


0.08 


0.025 


[Nev] 14.32 pm 


0.071 


0.20 


0.18 


0.10 


0.14 


0.012(7) 


0.009 


0.012 


0.012 


0.009 


0.014 


[Nem] 15.56 yum 


0.29 


2.2 


1.2 


3.0 


0.47 


0.05(3) 


0.009 


0.03 


0.03 


0.05 


0.006 


[Fen] 17.94 pm 


0.16 


0.90 


1.3 


1.8 


0.87 


0.03(2) 


0.005 


0.015 


0.03 


0.04 


0.013 


[Sm] 18.71 yum 


0.82 


5.0 


5.2 


2.5 


4.2 


0.08(2) 


0.025 


0.08 


0.15 


0.05 


0.06 


[Oiv] 25.91 yum 


2.3 


13 


11 


12 


7.2 


0.20(5) 


0.08 


0.24 


0.26 


0.26 


0.13 


[Sm] 33.48 yum 


0.28 


1.9 


2.3 


0.76 


0.6 


0.04(1) 


0.011 


0.05 


0.06 


0.015 


0.015 


[Sin] 34.81 yum 


1.0 


3.7 


3.9 


1.7 


1.7 


0.07(3) 


0.04 


0.10 


0.11 


0.03 


0.04 



Notes. The table contains data for the strongest lines only. The regions are numbered as in Hurford & FesetJ d 19961) with their positions shown in 
Fig. [7] The "average" column assumes single optical line reddening correction of 200 times. Note that the pixel sizes are different for the IRS LL 
(5'.'08 or 6.06 x 10" 10 sr, A > 15 yum) and SL (1'.'85 or 8.04 x 10~ H sr, A < 15 yum) modules. 



Th e rO ml far-infrared line r atio 7(51.81 yum)/7(88.36 pm) 
in the Borkow ski & Shuffl (Qj90) models is significantly higher 
than observed (6 versus 2.5). This discrepancy arises because 
they do not include emission from the pre-shock photoionized 
region. In the SD-200 model, in contrast, the pre-shock re- 
gion dominates the [O in] infrared line emission and diminishes 
this infrared line flux ratio to 1.2. Therefore, to reproduce the 
observed line ratio the pre-shock and post-shock contributions 
should be comparable (see also Table|6]i. 

The model I-H provides the most accurate description of the 
neutral oxygen FIR line ratio. In this model, most of the [Oi] 
emission originates in the post-shock PIR. The line ratio corre- 
sponds to the high-density limit, but the predicted line intensi- 
ties are about 20 times stronger than observed. As mentioned in 
Sect. 14.1.21 this discrepancy may be reduced if we assume higher 
pre-shock atom number density than in the model (for example, 
300 crrr 3 instead of 30 cm -3 ), as the line emissivities decrease, 
but both the ionization parameter and the column density of the 
PIR remain the same. 

3.4. Energetics of the infrared emission 

From the analysis of the Spitzer archival observations, it is 
straightforward to determine that in the wavelength range from 
5 to 35 pm the total Cas A luminosity is 3.3 x 10 37 erg/s (the 
luminosity values are compu ted assuming a distance to Cas A 
of 3.4 kpc; lReed et a l. ( 1995)). The total spectral line luminosity 
in the same range is found to be one order of magnitude lower, 
contributing 3.0 x 10 36 erg/s. 

For the longer wavelengths between 35 pm and 1 mm, we 
estimated the tota l Cas A luminosity using data published by 
iHines et all (120041) to be about 1.5 x 10 37 erg/s. The line con- 
tribution to this value in the narrower ISO LWS range between 
40 and 200 yum is about 9 x 10 35 erg/s, as follows from Table|3] 
Therefore, the total Cas A dust continuum and line luminosities 
in the infrared range between 5 pm and 1 mm are 4.7 x 10 37 and 
3.9 x 10 36 erg/s, respectively. 

It is interesting to compare these values with the X-ray and 
radio luminosities. In radio, the total luminosity at wavelengths 
longer than 1 mm is about 1.0 x 10 35 erg/s as ca l culate d di- 
rectly from the sp ectrum provided by iBaars et al] d 19771) and 
IHines et alJ ((2004), i.e., much lower than in the infrared. 



It is not so straightforward to determine the total Cas A lu- 
minosity in X-rays because of interstellar and internal X-ray ab- 
sorption. Various estimates res ult in values corrected for the ab- 
sorpti on of (1 -5) x 10 37 erg/s (Hartma nn et all[l997l:IZombeckl 
2007), but are dependent on the assumed spectral model and ab- 
sorbing material column density. 

It can be seen from this comparison that the infrared contin- 
uum emission is as an important radiative energy loss mecha- 
nism from the supernova remnant as the X-rays. 

It is easy to notice that the infrared line-to-continuum emis- 
sion ratio in the Spitzer range is variable over the supernova 
remna nt, as indicated also by the spectra provided by [Rho et alJ 
(2008), with a typical value of 0.10. The highest values of this 
ratio of about 0.22 are characteristic of the central Cas A region. 
In these regions, the energy losses in the line emission are rela- 
tively twice as high as average, although, still significantly below 
the losses in the infrared continuum. 

4. Physical conditions and abundances in the FMKs 

The measured line fluxes may be used directly to determine the 
ionic abundance^ if the model of the emitting region (temper- 
ature and density distribution) is known. Sicne the models do 
not reproduce observations precisely enough (see Table [T]i, we 
cannot assume that this distribution is known for the FMKs. 

Thus, the line fluxes cannot be used in such a direct fashion 
and the inferred ionic abundances will depend on the assumed 
physical conditions. Fortunately, emissivities of some of the fine- 
structure lines are only weakly dependent on the electron densi- 
ties and temperatures accross a rather wide range, allowing us 
to use these lines for the abundance determination even without 
detailed knowledge of the line-emitting region properties. 

An analysis of the measured flux ratios of lines of the same 
ion is a more powerful method. These line ratios may be used 
directly to constrain the temperature and density of the emitting 
region, provided that the emitting region is uniform and physi- 
cally the same for both lines defining the ratio. 

However, these assumptions are not fulfilled in the fast- 
moving knots of most line pairs and we cannot obtain emitting 

9 In this paper, by abundance we denote a ratio of number density 
of given ions or atoms to the number density of all oxygen ions (i.e., 
n/n ). 
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region parameters from, e.g., a comparison of optical and in- 
frared lines of [Om]. Even the density-sensitive ratios of two 
infrared lines should be interpreted carefully, as from the theo- 
retical models it follows that their emission may originate both 
in the pre-shock and post-shock regions that have comparable 
emission measures of the order of 10 19 cm" 5 and temperatures 
of the order of (0.5 - 2) x 10 4 K, but densities that differ by a 
factor of a thousand or more. 

More specifically, the pre-shock photoionized region has a 
low electron density (of the order of 100 - 250 cirT 3 according 
to the SD-200 model) and produces lines with approximately 
low-density limit line ratio, and the post-shock cooling region 
(n e « 10 6 cm -3 ) produces lines having high-density limit line 
ratio (i.e., the corresponding level populations are determined 
by collisional processes). Therefore, instead of determining the 
density, some of the FIR line ratios determine the contributions 
of the high- and low-density regions to the total line emission. 

In the following subsections, we analyze the line flux ratios 
to infer the physical conditions and the line fluxes to estimate the 
ionic and elemental abundances. 



4. 1. Information obtained from the same ion line flux ratios 

In the infrared range probed by the Spitzer and the ISO obser- 
vatories, there are seven line pairs that can be used for plasma 
diagnostics, and have at least one of lines defining the pair de- 
tected. These line pairs correspond to transitions between three 
of the lowest fine structure states 3 Po,i,2 in the ground p 1 and 
p 4 electron configurations. Of these seven diagnostic line pairs, 
three have both lines detected ([Om], [Sm], and [Arv] lines) 
and four have detections of just one of the lines ([Oi], [Nem], 
[Ne v], and [Arm] lines). The measured values of the flux ratios 
and the inferred physical parameters are summarized in Table [6] 



4.1.1. [O in] FIR line flux ratio 

The measured ratios of the far-infrared [Om] lines 
7(51.81 yum)//(88.36 pm) in different ISO spectra lie be- 
tween 0.7 and 2.5 (see Table [3]) with some indication of a 
correlation between the line ratio and the optical line intensity 
(a correlation is also observed for the [S m] line ratio map 
compiled from the Spitzer data). 

In the approximation of two emitting regions (pre-shock and 
post-shock), we use the following expression to determine the 
pre-shock PIR contribution to the total line b flux from the mea- 
sured a and b line flux ratio R = I a /h ( see Appendix iBl for its 
derivation): 



PIR,2> = 



'PIR.i 



R 2 -R 
Ri -R\ 



(4) 



where /pir^ and I lot ^ denote the photoionized region and total 
line b fluxes, respectively, and R\ and 7? 2 denote theoretical line 
ratio values in the pre- and post-shock regions. 

The [O m] line ratio of 2.5 observed for the ISO LWS region 
#2 with many optical knots may thus be achieved only if 85% 
of the 88.36 pm line and 40% of the 51.81 pm line originate in 
the pre-shock region. The corresponding theoretical line ratios 
in this case are Ri -1.1 andi?2 = 9.7 (se e Table [6). 

As shown by ISufherland & Dopital (1 19951) . the 5007 A 
line mostly originates in the pre-shock region. We determined 
that this is also the case for the 88.36 pm line. This al- 
lows us to take one further step and use their flux ratio of 
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Fig. 8. The constraints on the pre-shock region emitting the 
[Om] lines from the average 7(88.36 pm)/I(5007 A) line flux 
ratio (solid line). The dashed lines correspond to the constraints 
if the line flux ratio is changed by 20%. The contour denotes the 
temperature-density relation in the region of the SD-200 model 
most intensely emitting [Om] lines. 



7(88.36 pm)/I(5007 A) « 0.10 to estimate the pre-shock physi- 
cal conditions in the region emitting these [Om] lines. 

This line ratio is sensitive to both density and temperature 
and in a single-temperature emission region model corresponds 
to n e of several hundreds cirT 3 for temperatures (1.0 - 1.5) x 
10 4 K (see Table|6]i. These values are very similar to, but some- 
what different from, ones given by the SD-200 model for which 
T e < 1.1 x 10 4 K and n e « 200 cm -3 in the region that most 
intensely emits both [Om] lines. The resulting observational 
constraint on the pre-shock line emission region is indicated in 
Fig.|8] 



4.1.2. [Oi] FIR line flux ratio 

It is more difficult to use the [O i] line ratio for the quantitative 
analysis, because only one of the lines is detected. From ISO 
observations of region #2, we obtain the 3cr lower limit to the 
7(63. 19 /mi)/7(145.5 pm) fine ratio of about 30. 

Some information may still be obtained about this lower 
limit: it corresponds to either extremely low temperatures below 
100 K or high electron densities above 10 s cm -3 (see Table|6]l. 

According to the theoretical models, there is essentially no 
neutral oxygen in the post-shock cooling region because of slow 
recombination processes. The models also predict the rapid ion- 
ization of neutral oxygen in the pre-shock PIR with a very low 
line ratio 7(63.19 /im)/7(5007 A) much less than 0.01 (the ob- 
served value is 0.07, see Table©. The pre-shock region in which 
the [O i] line originates, is also excluded from the observed value 
of the 7(63.19 yum)/7(145.5 pm) line ratio. 

From this, we conclude that the far-infrared [O i] line emis- 
sion in the FMKs originates in the dense photoionized region 
after the shock. 

All the models predict values of this line flux ratio with re- 
spect to the 5007 A line to be much higher than observed (see 
Table [TJ. Part of this discrepancy may be explained by the low 
mo del values for the pre-shoc k oxygen densities. For example, 
the iBorkowski & Shull d 19901) models assume a very low pre- 
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Flux 


Single region 




Two regions 




OUCC LI LI 1 1 1 


line a 


line £ 




r c , K * 


n e , CffT 3 


R, 


R 2 


/PIR.O 


-f 

JPIRM 


ro t1 

IV l i 




14^ ^ um 
1H-J.J jUIIl 




< too 

r\ r\f\ 1 , , 1 n4 

200 - 1 x 10 


any 

> 2 x 10 5 










IV m J 


S 1 SI inn 
Jl.Ol jUIII 


oo.jo pm 


Z.J 


300 - 1.5 X 10 


C AA 1 AAA 

500 - 1000 


1 i 


9.7 


0.40 


U.OJ 


LUmJ 


So. Jo yum 


jUU/ A 


U.1U 


1.0 x 10 4 
1.5 x 10 4 


800 
200 






1.0* 
1.0* 


1 n* 
l.U 

1.0* 


TNp ml 


T.f. 01 urn 


1 5 ^(S um 


^ OR 
^ u.uo 


any 


any 










[Nev] 


24.32 pm 


14.32 pm 


< 0.2 


< 3 x 10 4 


> 3 x 10 4 










[Sm] 


33.48 /jm 


18.71 /im 


0.5 


600 - 3 x 10 4 


1500 - 3000 


1.8 


0.09 


0.85 


0.25 


[ Ar in] 


21.83 pm 


8.991 pm 


<0.10 


any 


any 










[Arv] 


7.914 pm 


13.07 pm 


1.6 


1500 - 3 x 10 4 


(4 - 10) x 10 4 


1.0 


6.7 


0.5 


0.9 


[Fen]"* 


1.3209 iim 


1.2946 pm 


0.75 


any 


> 2 x 10 5 










[Fen]"* 


1.3209 pm 


17.94 pm 


0.024 


1700-3000** 
1700-1800** 


any 

> 1 x 10 5 * 











* Assumed value 
** Derived value 

*** The [Fen] line flux ratios are discussed in Sect. 15. II 



shock oxygen density of 1 cirT 3 . Taking it higher by two orders 
of magnitude will decrease the post-shock 63.19 yum line inten- 
sity significantly. However, this is still not enough to reconcile 
the model predictions with the observations. The most probable 
reason for the low observed 7(63.19 pm)//(5007 A) line ratio is 
a combination of the post-shock region truncation and high pre- 
shock density, as discussed below in Sect. 15.51 

It should be noted that based on ISO data we cannot exclude 
part of the 63.19 //m line originating in a cold (T < 100 K or 
so) diffuse medium outside the optical knots (although, see con- 
siderations in Sect. 15.5b . This possibility may be fully assessed 
only by high angular resolution observations, which only new 
generation observatories such as Herschel are able to perform. 

4.1.3. The [New], [Arm], [Nev], [Arv] and [S in] line flux ratios 

All these lines are detected in the Spitzer data and, because of its 
much higher angular resolution, provide much more information 
for the analysis. 

The upper limits to the [Ne m] and [Ar m] line ratios (see 
Table[6]for details on these and other line ratios) do not constrain 
the plasma parameters much, since both o f them correspond t o a 
wide range of densities and temperatures dRubin et all 2001). 

Limits to the [Ne v] line ratio also do not provide strong con- 
straints part ly because of spe ctral fringing near 20 fim in the 
IRS spectra dEnnis et alll2006l) that increases the effective noise 
level. Nevertheless, the limit to the [Ne v] 24.32 pm line implies 
that the [Ne v] emission mostly or entirely originates in the post- 
shock cooling region. 

As for the case of [Om], we derived the fraction of the fine- 
structure line emission originating in the pre-shock region for 
[S in] and [Arv] lines. Most of the [S in] line at 33.48 /jm orig- 
inates in the pre-shock region, which we use to estimate the S 
abundance in Sect. 14.21 

Existence of Ar 4+ in the pre-shock region (see Table [6]l is 
understandable, given its relatively low Ar 3+ ionization potential 
of 60 eV, comparable to the ionization potential of 2+ (55 eV). 
Therefore, similar amounts of 2+ and 3+ in the pre-shock pho- 
toionized region in the SD-200 model give us a reason to deduce 
that Ar 3 + and Ar 4+ ions also ha ve similar number densities. We 
note that lChevalier & Kirshnerl(ll979l) demonstrated that the vis- 



ible [Anv] lines originate in the low-density environment, which 
we associate with a pre-shock PIR. 

4.2. Abundances from the flux ratios to the 5007 A line 

In the general case, ionic and elemental abundances derived from 
the individual line intensities are strongly dependent on the the- 
oretical model of the fast-moving knots. This model dependence 
may be minimized if the spectral lines are known to be emitted 
in the pre-shock or the post-shock region only. 

Some dependence on the underlying model still remains, but 
the model features used in our analysis are relatively robust. For 
example, the peak of the post-shock plasma emission measure 
around 10 4 K is determined by the cooling rate that slows down 
at 2 x 10 4 K as resonance ultraviolet lines of various ions stop 
being effectively excited. The low-temperature limit of the line- 
emitting region is computed from the upper transition level ex- 
citation energy, which is a simple but reliable estimate. 

The ionic abundances in the pre-shock region are stronger 
dependent on the underlying model, since we estimate them for 
all elements from average ionization potential of the oxygen ions 
in the SD-200 model. Thus, the pre-shock abundances may be 
incorrect by up to a factor of a several. 

As follows from the theoretical models, spectral lines of 
singly ionized species (e.g., [Arn], [Sin], etc.) arise mostly 
in the photoionized regions, whereas ones of highly-ionized 
species (e.g., [Ne v], [Mg v], etc.) arise mostly in the post-shock 
cooling regions. The only exception is the [Arv] line, but the 
pre-shock region contribution to its lines have already been esti- 
mated above. 

Below we discuss these two line groups separately and sum- 
marize the results in Tables[7]and[8] 

Because of extreme chemical inhomogeneities in the Cas A 
supernova remnant, one should not expect abundances in X-ray 
emitting low-density ejecta to be the same as in much denser 
optically-emitting knots. In our studies of infrared lines from 
dense knots, the abundances determined from optical spectra 
should be more reliable. However, optical abundance precision is 
not expected to be superior to a factor of two, because of the very 
approximate physical models used in the data analysis. Optical 
data are als o not available for Si and are especially uncertain for 
Fe and Mg (IChevalier & Kirshnerl fl978h . 
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Table 7. Estimated FMK post-shock cooling region ionic abundances «i on /"o- See text for details. 



Line 


///(5007 A) 


T C ,K 


e, cm 3 /s 


At, 10 3 s 


nim/fk) 


«ion/«o (model)* 


[Six] 1.430 yum 


0.0016 


2 x 10" 


2.9 x 10~ 8 


40 


7 x 10 -5 


5 x 10 4 


[Si vi] 1.964 fan 


0.004 


1 x 10 4 


3.1 x IQ- 9 


40 


0.002 


0.01 


[Mgv] 5.608 yum 


< 0.003 


6 x 10 3 


9.0 x 10-" 


30 


<0.002 


4 x 10~ 4 


[Arv] 7.914 /im 


0.005" 


1 x 10 4 


1.3 x 10~ 8 


60 


8 x 10~ 4 


5 x 10 s 


[Nev] 14.32 /im 


0.015 


1 x 10 4 


2.4 x IQ- 9 


50 


0.06 


0.002 



* Estimated based on SD-200 model. The elemental abundances assumed to deduce the model ionic abundances with respect to oxygen are given 
in Tabled 

** This line is partly arising in the pre-shock region. As follows from Table[6] the post-shock region contribution is ///(5007 A) ~ 0.0025. 

Table 8. Estimated FMK pre-shock ionic abundances ni on /"o computed assuming that the lines originate only in the pre-shock PIR. 
See text for details. 



Line 


///(5007 A) 


T e ,K 


n e , cm 3 


e, cm 3 /s 


IP, eV 


Af, s 


«ion/»(0) 


"dement/MO) ** 


[Am] 6.985 yum 


0.5 


2.0 x 10 4 


100 


4 x 10~ 8 


27.6 


5 x 10 b 


1.0* 


0.005 


[Arv] 13.07 fim 


0.003 


7000 


250 


8 x 10~ 7 


75.1 


5 x 10 7 


8 x 10~ 5 


0.005 


[Nen] 12.81 yum 


0.10 


1.5 x 10 4 


150 


5 x 10- 8 


41.0 


1.4 x 10 7 


0.06 


0.02 


[Fen] 17.94 //m 


0.03 


2.0 x 10 4 


50 


1.1 x 10~ 8 


16.2 


9 x 10 5 


5.5* 


0.005 


[Sin] 33.48 /mi 


0.04 


1.5 x 10 4 


150 


3.0 x 10- 7 


34.8 


1 x 10 7 


0.015 


0.05 


[Sin] 34.81 /an 


0.07 


2.0 x 10 4 


50 


1.5 x 10~ 7 


16.4 


1 x 10 s 


0.7* 


0.05 


[Om] 88.36 yum 


0.10 


1.0 x 10 4 


200 


4.0 x 10~ 8 


54.9 


3x 10 7 


1.1 


1.0 



* Shown to arises mostly in the post-shock photoionized region; see Sect. [5] 
** Atomic number densities with respect to oxygen, see Table [9] 



For our analysis, we summarized the available data on abun- 
dances in the Cas A ejecta (see Table |9]l and also used some 
"average" values. However, as noted above and clearly seen also 
from Table [9] the observed abundances vary dramatically from 
place to place. 

4.2.1. High-ionization ionic species from the post-shock 
cooling region 

An easy estimate of the plasma composition may be made if the 
problem is simplified by assuming that emission originates in a 
single layer with given temperature and electron density. In this 
case one can estimate relative number densities of two types of 
ions from their line fluxes. The line flux ratio is then expressed 
from Eq. (f3J as 

^ _ Vi Si » ion ,i At! 

h v 2 s 2 n ioni2 Af 2 ' 

where n; on is ionic number density, Af is the line emission dura- 
tion in the post-shock cooling plasma element, and it is assumed 
that the electron densities are the same for regions emitting both 
lines. The line emissivities entering Eq. (f5J) are computed as de- 
scribed in Sect. l2.2l 

We note that the application of Eq. (|5]l is possible in the SD- 
200 model because both the ion and electron density remain ap- 
proximately constant at T 5; 10 4 K. This, combined with the 
relatively weak dependence of the fine-structure line emissivity 
on temperature, allows us to use this simple estimate without 
introducing errors of more than a factor of two or so. 

Equation (0 may be inverted to determine the relative ionic 
abundances from the observed line ratios. In it, we always as- 
sume the second line in the ratio to be the [O ra] 5007 A line. 

In the following, we use the SD-200 model with the 
post-shock 2+ abundance in the line emitting region of 
n(0 2+ )/n(0) ~ 1%. We also take into account that only about 



2% of the total 5007 A line emission originates in the post-shock 
plasma according to the SD-200 model. This fraction is uncer- 
tain by up to a factor of about two, thus introducing a systematic 
error of similar magnitude in all post-shock ionic abundance de- 
terminations. 

There are no published results about the theoretically pre- 
dicted abundances of ions other than oxygen in the post-shock 
cooling region of the FMKs. To estimate them, we utilized 
the electron and ion densit ies and the cooling function from 
Sutherland & Dopiti (119951) and traced the ion ionization and 
recombination processes in the post-shock cooling phase us- 
ing a tomic data compil ed by Dima VerneiFI for the Cloudy 
code dFerland etallll998h . 

Another source of uncertainty in the ionic abundance deter- 
mination using our method are the model-dependent values of 
At, which were also computed by tracking the plasma cooling 
and recombination. However, calculations show that for different 
ions and lines the Af values differ by no more than approximately 
a factor of two, as determined by the general emission measure 
behavior. The resulting model abundances are compared with the 
observational results in Table [7] 

It is seen that the abundances inferred from the SD-200 
model seem to be generally consistent with the results of obser- 
vations, showing that the model in general correctly represents 
the temperature and density evolution in the post-shock cooling 
region. The [Arv] and [Nev] lines are the only ones showing 
significant differences from the model expectations, although the 
derived ionic abundances do not significantly exceed the corre- 
sponding elemental abundances. 

We note that the two highest ionization fine-structure lines of 
[Si x] and [Si vi] observed from the fast-moving knots have not 
been o bserved by Spitzer , but detected in the near-infrared spec- 
tra by Gerard v & Fesenl d2001l) . We determined their line flux 



http : //www . pa . uky . edu/~verner/atom . html 
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Table 9. Summary of elemental abundance n e iement/"o values in Cas A ejecta from various sources. The "average" data used in this 
paper are shown in the last column. 





CK78 


CK79 


Wil02 


LH03 




Laz06 




Used in this paper 


Method 


Optical 


Optical 


X-ray 


X-ray 




X-ray 






Ne/O 


0.025 


<0.013 


0.02 










0.02 


Mg/O 


<0.04 


<0.005 


0.004 






0.005- 


•0.03 


0.005 


Si/O 






0.025 


0.003- 


-0.015 


0.008- 


•0.06 


0.05 


S/O 


0.06 


0.002-0.7 


0.014 






0.006- 


•0.05 


0.05 


Ar/O 


0.004 


0.003-0.1 


0.0035 










0.005 


Fe/O 


0.003 


<0.001 


0.006-0.02 


0.002- 


-0.01 






0.005 



Data source abbreviatio ns: IChevalier & Kirshnej j!978h - CK78, IChevalier & Kirshnei] j!979h - CK79, IWillingale etail j2002h - WU02, 
Lamina & Hwani d2003h - LH03. lLazendic et al.1 faOOfli) - Laz06. 



ratios to the [O in] 5007 A hne using reddening measurements 
from the same regions bv lHurford & Fesed ll 19961) . 



4.2.2. Low-ionization ionic species from the pre-shock PIR 

The post-shock plasma recombination is significantly slower 
than the cooling and, even at T e » 300 K, the post-shock cooling 
plasma still mostly contains multiply-charged ions, resulting in 
essentially no infrared line emission from singly ionized atoms 
at wavelengths shorter than several tens of pm. 

In this section, we neglect any contribution to the line fluxes 
from the post-shock photoionized region. Its possible existence 
and physical conditions are discussed in Sect. [5] 

We can then write Eq. © in a form similar to Eq. (|5j for the 
pre-shock region given by 



h _ Vl S\ » e .l »ion,l Af] 
h V 2 £ 2 «e,2 «ion,2 Af 2 ' 



(6) 



The equation also contains the electron density ratio, because 
ions of different ionization potentials exist in regions of different 
electron densities. 

In contrast to the post-shock case, here the timescales At are 
determined by the ionization rates and may vary by an order of 
magnitude or more from ion to ion. 

From the SD-200 model, we adopt the average ionization 
potential as a function of time and compute the ionization state 
distribution of all elements with the simple approximation that 
this average ionization potential is equal for all of them. We also 
take the temperature profile from the model, which has some ef- 
fect on the resulting emissivities. Fortunately, precise tempera- 
ture values have only a minor effect on the obtained abundances, 
since the infrared line emissivities depend weakly on T e above 
10 3 K or so. 

In the pre-shock ionization front of the SD-200 model, each 
initially neutral atom is consequently ionized with maximum 
ionic abundance of each ion of the order of 50-100%. Therefore, 
the ionic abundances obtained in Table [8] should be quite close 
to the elemental abundances, provided that the At values are es- 
timated correctly and the pre-shock PIR provides the major con- 
tribution to the line intensitvFT. 

As a test case, we also included in Table[8]the similarly com- 
puted 2+ abundance from the [O in] 88.36 pm line. The result is 
very close to the expected value of one (even too close, given our 



11 Abundance of Ar 4+ estimated from [Arv] lines is much lower than 
the Ar elemental abundance, presumably because only some fraction of 
Ar atoms is ionized as far as Ar 4+ in the pre-shock region. 
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Fig. 9. Energy diagram of the lowest levels of Fe + ion with tran- 
sitions used in our analysis and corresponding line wavelengths. 
The electron configurations are given below the diagram; num- 
bers besides the levels deno te total angular momentum /. Data 
from Ikalchenk o et all d2007l) . 



uncertainties of about a factor of two), illustrating the reliability 
of our method. 

Although there are several [Fe n] spectral lines in the infrared 
spectral range probed by Spitzer that originate in the low-lying 
electronic configurations, the two most promising are blended 
with other stronger lines (one at 25.99 pm is blended with the 
[O iv] line at 25.91 pm and one at 35.35 pm is blended with [Si n] 
line at 34.81 pm). Therefore in the Spitzer wavelength range, we 
used for our analysis only the [Fe n] line at 17.94 pm (see Fig. [9] 
for the energy levels of Fe + corresponding to these transitions). 

From Table [8] it is seen that the Ar + , Si + , and Fe + abun- 
dances estimated assuming that their fine-structure lines arise 
exclusively in the pre-shock region are much higher than the cor- 
responding elem ental abundance values obta i ned from optical 
and X -ray data (IChevalier & Kirshneri Il979t IWillingale et all 
2002). The possible solutions of this discrepancy are discussed 
in Sect. O 
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5. Post-shock photoionized region 

5.1. On the too high pre-shock intensities of some lines 

Abundances of Ne and S derived assuming that the fine-structure 
lines of [Ne n] and [S m] are emitted in the pre-shock photoion- 
ized region of the SD-200 model agree reasonably well with the 
results of optical and X-ray observations. Simultaneously, simi- 
larly derived abundances of other singly-charged ions (Si + , Ar + , 
Fe + ) are unreasonably high. They are even much higher than 
their estimated elemental abundances (compare with Tabled- 

Within the frame of the SD-200 model, it is difficult to in- 
terpret these results and at the same time to comment on the 
absence of some other infrared lines, such as [Fern] line at 

22.93 fim and many others. 

There are several possible solutions to this mismatch be- 
tween the model predictions and observations. 

- The inconsistency might result from the application of the 
SD-200 model, describing the plasma dominated by C, O, 
Ne, and Mg, to the plasma, containing considerable amounts 
of S, Si, and Ar. However, the total abundance of these ele- 
ments with respect to oxygen reaches only about 10%, which 
should not significantly alter the shock structure. 

- Flux in the Si and Fe infrared lines might mostly originate 
in physically separated Fe-dominated clouds. This hypothe- 
sis is based on two observations: (a) in [Fen] line observa- 
tions Cas A has a differen t morphology than in the [O iv] line 
( Ennisetal l l200aiRholtan, l2008): and (b) X-ray Chandra 
observations indicate that hot Si- and Fe-dominated blobs 
exist in Cas A. 

However, the Cas A morphology in the Fe Ka X-ray line 
also has a strongly different morphology from that in the 
[Fen] 17.94 jjm line weakening these arguments. This hy- 
pothesis could be tested by observations of the Cas A ejecta 
with higher spectral resolution, which would allow one to 
determine whether [Fe n] and other infrared lines from the 
same knots have the same peculiar velocities. 

- The [Sin], [Arn], and [Fen] lines might originate in sep- 
arate weakly-ionized photoionized region of high emission 
measure. It might be situa ted e i ther be fore the shock (the 
"far precursor" region of Etoh| dl981alblV) o r after it (the 
post-shock PIR present in models of Itohl dl981albh and 
iBorkowski & ShuTO dl990h ). This hypothesis is discussed 
further below. 

It is possible to distinguish between lines emitted by the 
pre-shock and post-shock PIRs because of their strongly differ- 
ent electron densities. For this purpose, we use additional con- 
straints on this line-emit ting region that can be obtained from the 
lGerardv&Fesenl(l2001h observations of the near-infrared [Fe n] 
lines (see Tabled. 

The [Fen] line ratio 7(1.3209 //m)//( 1.2946 //m) < 1 (see 
Fig.|9]for the transition diagram between the Fe + levels) strongly 
indicates the electron density above 10 5 crrT 3 at any tempera- 
ture, thus weakening the hypothesis of the "far precursor". In 
the case of lower densities, this ratio is always higher, reaching 
values of 10 - 20 in the low-density limit because of strong dif- 
ferenc es in the excitation rates o f the respective transition upper 
levels dRamsbottom et all 120071) . 

Another problem for the pre-shock "far precursor" hypoth- 
esis is the observed ratio of the [Sin] 34.81 fim to the [Fen] 

17.94 /im lines of about 2. Both neutral and singly ionized Si and 
Fe atoms have similar ionization potentials, therefore we might 
expect that their ionization fractions are similar. From the line 



emissivities for the low-density conditions and average abun- 
dance ratio of n(Si)/n(Fe) « 10 (Table [9j, the expected line ratio 
is as high as 100. We therefore conclude that the PIR contribut- 
ing to the emission in these lines is likely to extend after the 
shock. 



5.2. Lines originating in the post-shock PIR 

All spectral lines of the neutral, singly-, and doubly-charged ions 
originate in the pre-shock and the post-shock photoionized re- 
gions. The strongly different densities of these regions (by a 
factor of more than 10 3 ) produces a much lower ionization in 
the denser post-shock PI R, which is only partially ionized (e.g., 
IBorkowski & ShullLll990h . 

In contrast, neutral and singly-ionized atoms are not abun- 
dant in the pre-shock PIR, where they are rapidly ionized to 
higher degrees. Estimates based on the SD-200 model confirm 
that for these atoms the pre-shock PIR contribution to flux in 
most of their lines is below 10-30%. 

For example, the [S n] lines close to 6720 A are strongly 
damped in the high-density medium of the post-shock PIR and 
are obviously emitted before the shock. It provides the possibil- 
ity of evaluating directly the pre-shock contribution to the bright 
near-infrared [S n] lines near 10 300 A from the low-density line 
ratios, which turns out to be about 10%. 

For our qualitative analysis of the post-shock photoionized 
region, presented below, we use the lines listed in Table [10] 
There are also other lines that originate in the same region, but 
their analysis does not yield additional conclusions. 



5.3. Charge exchange processes in post-shock PIR 

Atomic and ionic spectral lines in a rarefied partially photoion- 
ized cold plasma are produced by 

- Collisional excitation from the ground state, 

- Radiative recombination onto excited states, and 

- Charge exchange reactions into excited states. 

The charge exchange process at low temperatures is partic- 
ularly effective in our case of partially ionized plasma, as the 
atom-ion charge exchange is not slowed by electrical repulsion 
of the reactants. Since the oxygen atom ionization potential is 
higher than those of C, Si, S, and Fe, the charge exchange reac- 
tions may neutralize oxygen simultaneously ionizing these ele- 
ments. 

It is interesting to note that the charge exchange pro- 
cess have never before been accounted for in the theoretic al 
FMK models. The pure-o xygen FMK models dltohL 1 1981 alibi: 
IBorkowski & ShullL 1 1 990t) do not even mention charge exchange 
between, e.g., neutral and doubly ionized o xygen as one of th e 
rel evant processes. Based on th e method of[Dopita et al. ( 1984), 
the I Sutherland & Dopital d 19951) computations do not include the 
charge exchange reactions because the plasma tempe rature is too 
low for them to occur, according to HasteJ d!962l) . However, 
their argument does not apply, since in our case there are pseudo- 
crossings of intermittent molecular ion potential energy curves. 

We base our estimates of the rates of the relevant charge 
exchange reactions on t he si mple procedure suggested by 
iPeauignot & Aldrovandil d 1986b accounting for considerations 
of IStancill d200ll) . Our estimates are presented in Table QT| 
with the spectral lines produced by a charge exchange into 
excited states (so-called cha r ge exchange exc itation, see, e.g., 



excited states (so-calle d cnar ge e xcnange excite 
Dalgarno ^Stembergl(fT982l) : lDalgarno1dl985l) '). 
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Spectrum 


X 


///(5007 A) 


Source 


Observed? 


Application 


rie I 


IUojU A 


n aha* 


OrUl 


ies : 


He abundance 


[Ci] 


8727 A 


0.002 


HF96 


Yes 


Post-shock PIR disruption 


[Oi] 


6300 A 


0.08 


HF96 


Yes 


C abundance 




63.19 yum 


0.07 


This paper 


Yes 


O ionization state, electron density, 












extent of the post-shock PIR 




145.5 yum 


<0.0024 


This paper 


No 


Electron density 


Oi 


7774 A 


0.003 


HF96 


Yes 


Extent of the post-shock PIR, O ionization state, 












reference for line ratios 


rSinl 


34.81 yum 


0.07 


This paper 


Yes 


Extent of the post-shock PIR 


[Si] 


25.25 fim 


<0.01 


This paper 


No 


S° abundance 


[Sn] 


10287-10370 A 


0.11 


HF96, GF01 


Yes 


Probable multi-r structure of the post-shock PIR 


[Arn] 


6.985 fim 


0.5 


This paper 


Yes 


Extent of the post-shock PIR, Ar ionization state 


[Fen] 


1.2946 fim 


0.001 


GF01 


Yes 


Electron density 




1.3209 yum 


0.0007 


GF01 


Yes 


Electron density and temperature 




17.94 yum 


0.03 


This paper 


Yes 


Electron temperature, extent of the post-shock PIR 




25.99 yum 


<0.15 


This paper 


No? 


Electron temperature 



* An upper limit to the line intensity intensity (blended lines) 



Note. References: HF96 - iHurford & Fesenl fl996). GF01 - iGerardv & Fesenl d2001h . 



Table 11. Charge exchange reaction rates relevant to the post-shock PIR in an oxygen-dominated medium. 



Reaction 


AE, eV 


6E, eV 


»cx, cm J /s 


Emitted lines; notes 


+ + C° ^ 0°('D) + C + 


2.36 


0.39 


1 x 10" 


[Oi] 6300, 6363 A 


+ + Si -> O + Si + (3s 3p 2 4 P) 


5.47 


0.15 


3 x 10" 


[Sin] 2329-2351 A 


o + + S° -> O + S + ( 2 P) 


3.26 


0.22 


3 x 10" 


[Sn] 4069, 4076, 6716, 6731, 10287-10370 A 


0+ + Fe° -> O + Fe + (excited) 


5.72 


many transitions 


1.7 x 10" 


multiple lines; rate from Rutherford & Vroom (1972) 



Notes. It is assumed that the atoms and ions are in their ground states (or excited fine-structure states) before the charge-exchange reaction. AE 
denotes the difference in the ionization potentials of atom and ion in their ground states, 6E denotes the same difference for outgoing atom and ion 
in specified states, and a C x ls the charge exchange rate estimate. Only excited state classification is shown. 



Sicne the charge exchange is by far the most rapid of colli- 
sional processes, the ionization state of various elements in the 
post-shock PIR is insensitive to the charge exchange rate values. 
Taking these rates lower by one order of magnitude produces 
only an insignificant change in the ionization states (see Fig.flOl 
for example). Roughly, the number of ionized oxygen atoms de- 
creases by the total number of C, Si, S, and Fe atoms, provided 

that «o > «C + «Si + »S + "Fe- 

The charge exchange reactions neutralizing singly charged 
ions of neon and argon in collisions with atoms of oxygen is 
estimated to be very s lowT^ (of the order of 1 x 10 -13 cm 3 /s, 
iLiu & Dalgarnol ([1996)), therefore the argon and neon ioniza- 
tion state in the post-shock PIR is controlled by ionization and 
recombination processes. Charge exchange of Ne and Ar ions 
with other neutral atoms (C, Si, S, Fe) is not important since 
these elements are mostly ionized (see above). This qualitatively 
explains the strong observed [Oi] line at 63.19 yum and the ab- 
sence of fine-structure [Sii], [S i], and [Fei] infrared lines. 

The weakness of the neutral oxygen recombination line at 
7774 A is also indirect proof of the charge exchange process 
operating in the post-shock PIR [3 The 7774 A line is tradi - 
tionally used to estimate the extent of this region (lltohl [1986), 
as the line intensity is proportional to the amount of oxygen re- 
combination (in the absence of charge exchange, also to ioniza- 



12 We are grateful to the anonymous referee for pointing out the source 
of the Ne and Ar charge exchange rates with oxygen. 

13 We are grateful to the anonymous referee for pointing out the use- 
fulness of the O i 7774 A spectral line. 



tion) events, which depends on the optical depth of the post- 
shock PIR. At the post-shock PIR temperatures, each radiative 
recombination of a singly ionized oxygen atom is followed by 
the 7774 A line emission with a probability of about p = 0.4 
dPequignot & Aldrovandlll986l) . 

However, since ionized oxygen in the post-shock PIR may 
be neutralized by both recombination and charge exchange, the 
7774 A line weakens, and its flux alone does not provide enough 
data to determine the optical depth of the post-shock PIR. To 
obtain full information about the number of recombinations in 
the post-shock PIR plasma, deeper observations are needed that 
would also detect optical and near-infrared recombination lines 
of other elements (e.g., Ci line at 10695 A, S i line at 9237 A, 
etc.). 

However, a lot of useful information may be extracted from 
the flux of the 7774 A recombination line. For example, one can 
roughly estimate the oxygen ionization from the ratio of its flux 
to the total flux from the post-shock PIR. 

If plasma contains only oxygen, then the number of O re- 
combination events will be equal to the number of absorbed ion- 
izing photons. The presence of other elements with similar ion- 
ization cross-sections and recombination rates will decrease the 
number of oxygen recombinations by a factor no/«t, if «t re- 
mains constant and no charge exchange reactions occur. In this 
case, the ratio of the 7774 A line flux to the total flux emitted 
from the post-shock PIR cannot be smaller than 

7(7774 A) = phvrmno = 0.4xl.6eV » o =0Q21 »o 
Z 1 /^ioniz «t 27 eV n t ' n t ' 
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0.5 1 1.5 2 2.5 

n(C)/n(e) 

Fig. 10. Ionization state of carbon in an oxygen-carbon par- 
tially photoionized plasma. Since carbon and oxygen photoion- 
ization and recombination rates are similar, the carbon ioniza- 
tion fraction is almost exclusively dependent on the ratio nc/n e . 
Dependencies on both total atomic density and carbon abun- 
dance at constant « e are minor. The upper curve corresponds to 
the charge exchange rate given in Table [TT] and the middle and 
lower curves to the rates artificially diminished by factors of 10 
and 100. It is seen that the charge exchange processes result in 
carbon being ionized by more than 95-99% when the number 
density of carbon atoms and ions is less than that of electrons. 

where /jVkmiz is the average ionizing photon energy. It is esti- 
mated to be 

_ 1 H sh _ 1 (5/2)k B T sh x(nJn t + l)n t /n 
2 N ph 2 N ?h 

where r s h and n e /n t are the temperature and the electron-to- 
atom ratio at the shock front, respectively, ss 2250 eV is the 
enthalpy per oxygen atom at the shock front (the cooling post- 
shock plasma should be isobaric most of the temperature range), 
and A^ph ~ 42 is the average number of emitted ionizing photons 
per oxygen a tom (its value being derived from the ionization pa- 
rameter Q of lSutherland & Dopital dl995lV ) in the SD-200 model. 
The factor of 1 /2 is the fraction of downstream photons. 

However, the observed value of this 7774 A line ratio is at 
least a factor of three lower (its ratio with respect to the FS line 
of [Am] alone is about 0.007, see Tables [TOl and 151). We note 
that this ratio of 1/3 is an absolute upper limit. The ratio of the 
expected 7774 A line flux to the observed one may be much 
higher, since part of the much stronger dust continuum emission 
(of flux 10-40 times higher than that in the 5007 A line) may 
originate in the post-shock PIR. 

From this, we infer that at least two thirds of the ionized oxy- 
gen atoms recombine by means of charge exchange. Assuming 
additionally that 

- this results in a complete ionization of atoms of other el- 
ements connected to oxygen by charge exchange reactions 
(see above), 

- the recombination rates of all singly ionized atoms are simi- 
lar, and 

- the average ionization de gree of plasma is about 40% 
(iBorkowski & Shul]Lll990h . 



we find that the total number of elements related to oxygen by 
charge exchange reactions constitutes about 2/3 of the number 
density of oxygen itself (i.e., «o/ n t ~ 0.6) and the oxygen ion- 
ization n(0 + )/«(0) is not higher than 0.10-0.15. 

5.4. The post-shock PIR temperature and density 

Once the physical situation has been established, we now deter- 
mine the conditions in the post-shock photoionized region. 

The lower limit to the temperature of about 1000 K is ob- 
tained from the requirement that the [Fe n] 25.99 fim line emit- 
ted by this region should not be too bright (an observational con- 
straint is 7(25.99 yum)//(17.94 /urn) < 5, see Sect. l3"XTT l. 

The very low value of the near-infrared [Fen] 1.3209 fim 
line ratio to the 17.94 fim line corresponds to T e « 1700 K in 
the high-density limit (the high electron density follows from 
the [Oi] fine-structure line ratio, see Sect. 14.1.21 . This [Fen] line 
ratio is very sensitive to the temperature because of the different 
excitation potentials of the transitions upper levels and possibly 
provides only an upper limit to the temperature with the near- 
infrared 1.3209 fim line partially arising in some smaller and 
hotter region. 

Assuming in our analysis the PIR temperature of 1500 K, 
we can obtain stronger constraints on the electron density from 
the limit on the FIR [Oi] line ratio, assuming that these lines 
originate in the post-shock PIR, as discussed in Sect. 15.51 

In this case, the FIR [Oi] line ratio corresponds to n e ~ 
5 x 10 5 cm 4 . We assume this value for the electron density 
and twice as large a value n t = 1 x 10 6 cm -3 for the total atom 
and ion numb er density (average ioniz ation of about 40-50% is 
obtained by IBorkowski & Shulll (Q990) and should be approxi- 
mately correct, since recombination rates of all singly charged 
ions are similar). 

We note that based on the assumption of constant post-shock 
gas pressure this n t value corresponds to the pre-shock num- 
ber density of about one hundred atoms per cm 3 , implying that 
the magnetic field does not significantly modify the post-shock 
plasma compression. 

5.5. The post-shock PIR ionization and thickness 

The ratio of the neutral oxygen lines 7(7774 A)//(63.19 fim) « 
0.045 may be used to infer the oxygen ionization degree, which 
is found to be «(O + )/«(O ) * 0.06 for our assumed electron 
density. This estimate is precise up to a factor of two if n e is in 
the range between 1 x 10 5 and 1 x 10 6 crrT 3 and only weakly 
depends on the region temperature. 

The value of the oxygen ionization degree is in reasonable 
accordance with the estimates of Sect. 15.31 This is also a reason 
for us to assume that most of the 63.19 fim emission originates 
in the FMKs, and not in some more diffuse regions around them 
(remember that this line was observed by ISO, which had poor 
angular resolution). In its turn, this allows one to use the [Oi] 
63.19 fim line intensity to estimate the extent of the post-shock 
PIR, assuming that it is homogeneous. 

To estimate the post-shock PIR thickness (see Fig. Q}, we 
can also use the fine-structure lines of [Sin], [Arn] and [Fen], 
as well as the O i 7774 A recombination line. These estimates 
are summarized in Table fT2l 

Closer agreement between the region sizes obtained from 
the [Oi] and [Fen] lines may be achieved if the region tem- 
perature is assumed to be about 1800-2000 K (then the [Fen] 
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Table 12. Estimates of the post-shock photoionized region thick- 
ness Zpspir from several fine-structure line ratios to the 5007 A 
line. 



Spectral line 


£, erg cmVs 


n(ion)/n(elem) 


^pspir, cm 


[Oi] 63.19 


1.2 x 1(T 24 


0.94 


1.2 x 10" 


[Sin] 34.81 fim 


1.3 x 1(T 23 


1.0 


2.1 x 10" 


[Arn] 6.985 yum 


2.0 x 1(T 21 


0.4 


2.0 x 10" 


[Fen] 17.94 //m 


4.0 x 10~ 23 


1.0 


3.0 x 10" 


O i 7774 A 


1.3 x 10~ 24 


0.06 


8.0 x 10 10 



Notes. Estimates made for post-shock PIR electron density n, = 5x 
10 5 cm" 3 and temperature T c = 1500 K. Oxygen abundance have been 
taken to be no/n t = 0.6 (see Sect. 15. 3b . other abundances as in Table[9] 



line emissivity will increase decreasing its region size estimate 
to (2.0- 2.2) x 10 11 cm). 

The estimate from the 7774 A line also becomes closer to the 
value derived from the [O i] line at these slightly higher temper- 
atures, becoming about 9.3 x 10 10 cm. Alternatively, decreasing 
the oxygen ionization degree increases the post-shock PIR size 
estimate obtained from this line. 

It is worth noting that values of the post-shock PIR size in 
Table Q~2] are only about factor of about two smaller than th e 
value of 3 x 10" cm predicted by Borko wski & Shulll d!990h . 
This observational evidence therefore does not require si gnif- 
icant truncation of the post-shock region, proposed by lltohl 
(1986). However, other spectroscopic signatures, discussed in 
Sect. 15. 91 show that some truncation of the post-shock PIR seems 
to be present. 

5.6. Amount of carbon in the FMKs 

The [O i] 6300 A line originates in neither t he pre- shock PIR, nor 
the po st-shock cooling region, as shown bv lSutherland & Dopital 
(119951) . therefore its source is located in the post-shock PIR. 
Since it is a factor of 25 stronger than the Oi 7774 A recom- 
bination line, we can conclude that recombination provides only 
a minor contribution to the line flux. It may be produced by col- 
lisional excitation, but then the corresponding temperature from 
the 7(6300 A)//(63. 19 yum) line ratio of about unity corresponds 
to temperatures of at least 3000 K. 

However, the 6300 A line may be produced by a charge ex- 
change reaction with carborQ ( see Table fTTt . As this reaction 
is very fast, almost all the carbon is ionized and every carbon 
recombination will be immediately followed by the charge ex- 
change reaction with an oxygen ion. Then, with a probability of 
about 72% (accounting for collisional transitions downwards), 
oxygen emits the 6300 A line. A rather strong observed 6300 A 
line infers a considerably high abundance of carbon in the FMK 
plasma. 

Specifically, from this line ratio to the 7774 A line, one can 
conclude that the carbon abundance is 

n c /«o ~ 0.5, 

assuming that all carbon is ionized. 

This result is only weakly dependent on the post-shock PIR 
temperature, as long as it stays below 3000 K, where collisional 
excitation of 6300 A begins to play a role, and the electron den- 
sity, as long as it is lower than about 5 x 10 6 crrT 3 , where col- 

14 Note that the charge exchange with other elements do not produce 
the 6300 A line; carbon is unique in this respect. 



lisional de-excitation in neutral oxygen starts to diminish the 
6300 A line emissivity. However, estimated carbon abundance 
is proportional to the assumed oxygen ionization, taken here to 
be 0.06 (see Sect. 15. 5b and is therefore uncertain by up to a factor 
of two. 

5. 7. Sulfur in the post-shock PIR 

From the non-detection of the [S i] fine structure line at 
25.25 /im, we can define an upper limit to the neutral sulfur abun- 
dance in the post-shock PIR. For this purpose, we compare the 
[Si] 25.25 //m line with the [Oi] 63.19 yum line having similar 
dependence on temperature and density. The upper limit on their 
ratio 7(25.25 //m)//(63.19 /mi) < 0.14 corresponds to the abun- 
dance ratio of «(S )/ n (O°) < 0.004, showing that sulfur is more 
than 90% ionized. 

We note that in the framework of our model it is difficult 
to understand the brightness of the [S n] lines near 10 300 A. As 
discussed in Sect. 15. 21 they are emitted after the shock. However, 
in our simple single-temperature model none of the line pro- 
duction mechanisms - radiative recombination, charge exchange 
with oxygen (see Table fTTI) . or collisional excitation - can ex- 
plain the high brightness of these lines. 

5.8. Helium in the optical FMKs 

X-ray data infers that a large a mount of elements w ith Z < 8 are 
present in the hot Cas A ejecta (Hughe s et al.Ll200C ). Obviously, 
these elements are hydrogen, helium or carbon. iDewev et aU 
(2007) even suggested that hot X-ray emitting plasma consists 
of mor e than 95% He. This i s also consistent with the sugges- 
tion by I Sutherland & Dopital (Il995h that a hot plasma between 
the optical FMKs would be much brighter in X-rays than ob- 
served, if it consisted mostly of oxygen. 

The amount of helium in the FMK plasma may be con- 
strained, e.g., from the upper limit to its well-known near- 
infrared recombination line at 10 830 A. Each singly charged he- 
lium ion recombination with a probability of about 14% results 
in the emission of this line (Peauignot et all fl 99 lh . Therefore, 
the total recombination of helium after the shock will result in 
0.14n(He)/n(O) photons per oxygen atom. 

Comparing this with the SD-200 model 5007 A line in- 
tensity of about 7 photons per oxygen atom, we would ex- 
pect the 7(10830 A)//(5007 A) * 0.02n(He)/«(O). The ob- 
served line ratio upper limit of 0.006 thus corresponds to the 
helium abundance upper limi t of about »(He)/»(Q) ^ 0.3 . This 
result is consistent with the IChevalier & Kirshnerl (fl978h and 
iPeimbert & van den Berghl (1 19711) estimates, but is 5 - 10 times 
more restrictive because of more certain theoretical model of the 
post-shock PIR. 

This value is probably underestimated by up to a factor of 
two as helium will still be partially ionized in the post-shock 
PIR and may not recombine fully if this region is disrupted (see 
below). Nevertheless, it is obvious that the helium abundance of 
optical FMKs is not very high, implying that optical and X-ray 
emitting ejecta indeed have different chemical abundances. 

5.9. Evidence of post-shock PIR disruption from spectral 
lines 

Post-shock region disruption have been observed in both numer- 
ical simulations of the shock interact ion with the dense inter- 
stellar clouds in supernova remnants (Patn aude & FeseriL 12005) 
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and in laser experiments emulating this interaction (IKlein et all 
I2003h . 

As shown bv litohl dl9S6h . there are spectroscopic signatures 
that allow one to determine if the post-shock PIR is disrupted. 
Above, we have argued that the original diagnostics of lltohl 
(1986) - the O i 7774 A line - is not applicable because of the 
charge exchange between oxygen and other elements, and the 
poorly constrained typical chemical composition of the FMKs. 

Instead we consider other spectral lines that are sensitive to 
the final recombination in the post-shock PIR. Indeed, if this re- 
gion is present, then each ion in the post-shock PIR recombines 
at least once, giving a lower limit to the line intensities due to 
recombination. The most sensitive to the presence of final re- 
combination are high-excitation lines of atoms that should be 
completely ionized due to charge exchange with oxygen, i.e., 
ones that should not be emitted in any other region of a FMK. 

As an example, we consider the observed [Ci] line at 
8727 A. It is emitted with probability of roughly 10%, when a 
carbon ion recombines. Thus, there should be at least 0.05 pho- 
tons in this line per every oxygen atom (using the carbon abun- 
dance from Sect. I5.6l l. Comparing this with the SD-200 model 
5007 A line intensity of about 7 photons per oxygen atom, we 
would expect to measure 7(8727 A)//(5007 A) * 0.007. The ob- 
served 8727 A line is a factor of three less intense, implying that 
not all carbon recombines after the shock. 

Unfortunately, this result, which confirms the post-shock 
PIR disruption, is not very robust, since the carbon abundance 
value is not very tightly constrained. However, similar diagnos- 
tics may be performed with this and other lines using future more 
sensitive observations. These observations will show us the max- 
imum extent to which each atomic species recombines in the 
FMKs after the reverse shock. 

We note that another argument exists for the post-shock PIR 
truncation, which is based on energy considerationSThe en- 
ergy of the post-shock PIR line emission per oxygen atom is a 
factor 5-10 less than the energy that should be deposited if the 
emission were optically thick to the post-shock cooling plasma 
emission. The most obvious solution to this problem is post- 
shock PIR truncation. 

However, another possibility is that some part of the cool- 
ing plasma energy is transferred to dust heating inside the post- 
shock cooling region or post-shock PIR and is then emitted in the 
infrared dust continuum, which in total emits a factor of 12 more 
energy than the lines from the post-shock PIR (see Sect. 13.41 ). 
Most probably, a combination of these two factors (disruption 
and continuum emission) explains the apparent difference be- 
tween the expected and observed post-shock PIR emission-line 
fluxes. 

5.10. Effects of the dust on the post-shock PIR structure 

Using the data from iRho et al ] (I2008I) . it is straightforward to 
estimate the dust mass in each pixel of the Spitzer map of the 
Cas A. Its comparison with the optical image of this supernova 
remnant allows us to determine the dust-to-gas mass ratio in the 
bright knots from the [O m] line flux predicted by the SD-200 
model. 

If the emitting dust with temperature 7^ w 100 K dRho et all 

2008) is distributed over the entire volume of the knot, then the 
dust-to-gas mass ratio is of the order of 2-5%. If, on the other 



This argument was pointed out by the anonymous referee. 



hand, the dust emits only in the line-emitting regions, then the 
dust-to-gas mass ratio is of the order of unity. 

If the latter is true, then the current shock models do not 
apply and construction of an entirely different model is required 
to describe the shock structure of the Cas A knots. This analysis 
is outside the scope of this paper. 



6. Recombination lines in the infrared range 

The metal recombination lines (RLs) are another good tracer of 
the cold ionized plasma. Their emissivities as functions of tem- 
perature have no exponential cutoff at low T e and increase with 
decreasing temperature approximately as T~ l at T e < 10 4 K. The 
line emissivities are only weakly dependent on density. 

However, their emissivities are 4-6 orders of magni- 
tude less than those of the fine-structure infrared lines at the 
low-density limit. The infrared RL emissivity ratios to the 
fine-structure lines increase with electron density due to colli- 
sional effects that diminish the fine-structure line emissivities 
as n~ l starting from some critical density. At post-shock elec- 
tron densities of 10 6 cm" 3 , th e emissivity ratios are much higher 
dDocenko & Sunvaevl 12008b . but the infrared RLs are neverthe- 
less still far less intense than the fine-structure lines. 

As two examples, we consider infrared RLs near 10 and 
60 pm. The brightest regions emitting RLs are the post-shock 
cooling and photoionized regions (due to their high emission 
measures n 2 e l) and cold region between the photoionization front 
and the shock (due to its extremely low temperature and large 
extent). 

The Oi 5a recombination line at 7.45 pm originates in 
the post-shock photoionized region and its strongest component 
emissivity of 5 x 10~ 27 erg cm 3 /s is about 200 times less than that 
of the O i 7774 A line. The expected recombination line flux is 
therefore 3 x 10 4 times less than that of the [Am] 6.985 pm line, 
or about 0.5% of the background continuum emission. However, 
other elements will also produce recombination lines at the same 
wavelength (e.g., ionized carbon is almost an order of magnitude 
more abundant than ionized oxygen, see Sect. 15.6b and the total 
5a line flux may be a factor of 10-20 more. Estimates show that 
this limit is achievable with Spitzer with exposure times of a few 
hours. 

Another example is a Oi 11a RL at 69.03 pm. Its emissiv- 
ity is approximately 5000 times less than that of nearby [Oi] 
fine-structure line at 63.19 pm. Expected intensity, if summing 
with 11a lines of other ions, is about 5 x 10~ 15 erg/cm 2 /s from 
the ISO region #2. This is approximately equal to the expected 
5cr 1-hou r point source sensiti vity of the Herschel PACS in- 
strument (Pogl itsch et all l2006h . Since the Cas A is not a point 
source at the Herschel resolution (PACS has a 10" pixel size and 
most of the Cas A emission is contained within 4-8 PACS pix- 
els), one needs at least a few hours to achieve the 5cr detection 
of this recombination line. 

The recombination lines of highly-charged oxygen ions from 
the post-shock cooling region are expected to be about an or- 
der of magnitude dimmer than the O i RLs because of the lower 
emission measure in the region emitting these lines. 

The metal recombination line observations in the mid- and 
far-infrared are generally more difficult than in the optical and 
near-infrared due to lower line emissivities and higher back- 
ground continuum emission. They are also unable to provide 
information about the ion producing the spectral line, but only 
about its electronic charge. However, these observations would 
allow us to determine abundances of ions residing in cold regions 
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that do not produce fine-structure lines. Special emphasis should 
be placed on the singly-charged ion lines, since their recombi- 
nation lines in the optical and near-infrared ranges are divided 
into many weak components, which makes them more difficult 
to detect dDocenko & SunvaevLkoOSl) . 

7. Conclusions 

We have performed an analysis of the supernova remnant Cas A 
fast-moving knot infrared line intensities, by comparing obser- 
vations with the predictions of various theoretical models, which 
describe the FMK emission originating from the reverse shock 
interaction with the pure-oxygen or oxygen-dominated clouds. 
For this comparison, we have analyzed archival observational 
data acquired by the ISO and Spitzer observatories. We con- 
clude that accounting for the electron conductivity is essential 
to reproducing the observed l ine ratios. This is the reason why 
the lBorkowski & Shulll ([l990) model BS-DC provides the most 
accurate description of the observed oxygen line relative fluxes, 
although it is only precise up to a factor of several (the [Oi] 
lines are overpredicted partially because of the too low model 
pre-shock density). 

This emphasizes the need for more sophisticated models de- 
scribing shock propagation in the oxygen-dominated plasma. In 
this article, many line flux ratios are derived that will help us 
to construct these future theoretical models. Future far-infrared 
and optical observations of higher sensitivity and angular reso- 
lution will allow one to obtain far more information about the 
fast-moving knots once the models are constructed. 

Analysis of the infrared lines of O, Ne, Si, S, Ar, and Fe ions 
have confirmed the existence of three regions contributing to the 
infrared and optical line emission: 

- Lines of the intermediate- and highly-charged ions that arise 
in the post-shock rapidly cooling region (e.g., [Ne v], [Arv], 
[Si vi], [Six] lines). 

- Singly-charged ions that radiate mostly from the post-shock 
photo-ionized region (e.g., [Oi], [Sin], [Am], [Fen] lines). 

- Weakly- and intermediate-charged ion lines originate in the 
pre-shock photoionized region as well. It gives major contri- 
bution to the emission of lines of [Ne n] , [O in] , [S m] , [Ar v] , 
etc. 

An analysis of the FIR observations jo i ntly w ith the 
iHurford & Fesenl d 19961) and iGerardv & Fesenl (1200 ll) optical 
and near-infrared data has allowed us to construct a consistent 
model of the post-shock photoionized region. This became pos- 
sible by accounting for the charge exchange process that dimin- 
ishes oxygen ionization in the post-shock PIR. Observed high 
(> 0.9) ionization of sulfur and low (« 0.06) ionization of oxy- 
gen confirms the importance of the charge exchange reactions in 
this region. 

From the observed line flux ratios, it then follows that the 
post-shock PIR has atemperature T e « (1500-2000) K, electron 
and ion densities of n e « 5 x 10 5 cirr 3 and « t « 1 x 10 6 cirT 3 and 
is approximately (1 - 2) x 10 11 cm thick, assuming that oxygen 
constitutes about 70% of all atoms in the FMK plasma. Given 
the uncertainties of the model, this result is consiste nt with the 
thickness predicted in the Borkowski & Shull (1 1990b model. 

We have constrained the helium abundance to be 
n(He)/«(0) < (0.3 - 0.6) and the carbon abundance to be 
about n(C)/n(0) ~ 0.5. 

However, we recall that the obtained values of elemental 
abundances and ionization fractions in the post-shock PIR are 



affected by unknown biases, as our assumptions (plane-shock, 
single-temperature model) clearly do not adequately represent 
the structure of this region. 

From the comparison of observed infrared line intensities 
with the theoretically expected values, we conclude that 

- The theoretical models correctly describe the general struc- 
ture of the FMKs; 

- Accounting for the electron conductivity brings the BS-DC 
model into far closer agreement with the observed oxygen 
FIR line ratios to the [On] 5007 A line than the BS-F and 
SD-200 models; 

- Accounting for the emission from the photoionized region 
before the shock front is essential to reproducing the FIR 
[Om] line ratio; 

- The FIR [Oi] line inten sities in the Htoh| d!981bl) and 
iBorkowski &Shull ( fl990T) models are partially overesti- 
mated because of the too low pre-shock atom number densi- 
ties; 

- The pre-shock atom number density is at least 100 cirr 3 and 
may be as high as 300 cm 4 ; 

- The high density of the post-shock PIR suggests that there is 
no significant magnetic-field pressure in this region; 

- Th e pos t-shock photoionized region truncation, as suggested 
by |ltoh| ([1986), may be needed to explain the weakness of 
several lines produced by the final plasma recombination. 

The infrared metal recombination lines are shown to be de- 
tectable by the planned far-infrared instruments and useful for 
the derivation of the plasma properties. 

The high brightness of the [S n] lines close to 10300 A re- 
mains unexplained by our single-temperature model of the post- 
shock PIR. 
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Appendix A: Spitzer data cube PSF size estimate 

As described by ISmith et al.l (120071) . the mathematical transfor- 
mations performed by the CUBISM software change the point 
spread function (PSF) of the final data cube away from the initial 
PSF of the Spitzer IRS. This effect should be more pronounced 
in images with only a few pixels per PSF width, as constructed 
by the Spitzer IRS data cubes in the second order of each spec- 
trograph module. 

To determine the PSF changes induced by the data cube re- 
construction, we fit a two-dimensional Gaussian function to the 
shape of a point-like source 2MASS 23233176+5853204 situ- 
ated in the map of the SL2 module. This fit indicated that (here x 
is the coordinate along the individual slits and y - across them): 

- The source x centroid oscillates with wavelength with an am- 
plitude of about 0.3 pixels. The y centroid remains constant 
to within 0.05 pixels. 

- The PSFs in both the x and y directions are larger than the 
Spitzer IRS true PSF, computed, e.g., by the stinytim soft- 
ware. The full width at half maximum (FWHM) in the x and 
y directions of the spectral maps are greater than the IRS true 
PSF by 0.55 pixels and 0.9 pixels. 

The presence of th ese features is infer red by the CUBISM 
algorithm description in Smith et al. (2007). Although the values 
of the increase in PSF dimensions may seem insignificant at first, 
they are often comparable with the extent of the IRS true PSF 
(e.g., at 26 pm, the initial FWHM of the PSF is equal to the 1.2 
LL1 module pixel). 

Not accounting for the described increase in the PSF size be- 
cause of processing by the CUBISM software in our case would 
result in the Spitzer maps being seemingly more diffuse than e.g., 
optical maps smoothed to the IRS angular resolution. 

Although we have measured these effects only in the SL2 
module, we assume that the same PSF broadening in pixels oc- 
curs in other modules as well. This conclusion is qualitatively 
confirmed by visual comparison of optical maps smoothed to 
the corresponding resolution of the data cube spectral line maps. 

Appendix B: Derivation of Eq. © 

For the approximation of two homogeneous emitting regions 1 
and 2, the observed line a and b flux ratio R = I a /It, may be 
expressed as 

_ _ h _ Ia,\ + h,2 _ h,\ h,l + h,2 L h,l \ 

h h,\+h,i h,\h,\+h,i h,2\ h,\+h,2i 
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Denoting the line flux ratios arising in regions 1 and 2 as Ri = 
Ia,ilh,\ and R2 = I a ,ilh,2, and introducing the fraction of the 
total line b emission arising in the region 1 as 

, _ h,\ 

Jb,l = 



4,1 + h,i 



we obtain 



R = R l f b , l +R 2 (l-fb,i)- 



Expressing ft, 1 from this linear relation, we finally derive 
Eq. ©: 

R2-R 



